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TECHNICAL MEMORANDUM 


THE IMPACT RESPONSE OF CARBON/EPOXY LAMINATES 
(MSFC Center Director’s Discretionary Fund Final Report, Project No. 94-13) 


1. INTRODUCTION 


The term “composite materials” can be used to describe a host of different material forms used 
for a vast variety of final applications. The one common denominator for rill of these materials is that two 
or more distinct phases of material are combined in a manner such that the resulting “composite” 
properties are more desirable for the application at hand than any of the constituents alone. One of these 
phases is usually called the matrix, which acts to hold together the composite material. The other major 
phase is commonly called the rei nf orcement. Composites generally fall into three major categories defined 
by the type of reinforcement phase: particulate, short fiber, and continuous fiber. 

Particulate composites consist of particles (aspect ratio approximately less than 10) bound together 
by the matrix phase. The particles usually do not act as reinforcements for the matrix, but are commonly 
used as fillers (such as rocks in cement to make concrete) or to enhance properties other than strength 
(such as rubber particles in a brittle plastic to enhance the plastic’s toughness). 

Short fiber composites contain fibers with an aspect ratio of approximately 10:10,000. With 
the exception of most ceramic matrix composites (which use fibers for toughening), the fibers are the 
major load-carrying component of the composite material. This class of composites represents the most 
widely used and common composites, mainly due to what is commonly referred to as “fiberglass.” 
Strictly speaking, fiberglass is glass that has been melted and drawn into fiber form, although the term 
“fiberglass” usually means chopped glass fiber in a polyester or other inexpensive resin matrix. This 
material has been widely used since the 1950’s for components such as boat hulls, automobile bodies 
for sports cars, and recreational products such as surfboards, skis, and racket frames. 

The use of continuous fiber-reinforced composites is growing at a rapid rate since the cost 
of materials and processing is decreasing. These types of composites make the most efficient use 
of the hi gh stiffness and strength of the fibers used in them, and are usually what is being referred to when 
the term “advanced composites” is used. Some of the more common fibers used in this category are 
glass, carbon, and Kevlar™. Glass fiber-reinforced plastics (FRP’s) have good tensile and compressive 
strengths and modulii but are heavier than carbon FRP’s, which offer very high specific properties but are 
relatively weak in compression, limiting them to use in such tensile load-dominated structures as filament 
wound pressure vessels. 



High strength can be achieved by aligning the fibers in the direction(s) in which the structural 
load(s) will be applied. However, for laminated (layered) composites (which are almost always plates 
or shells), there are no fibers in the through-the-thickness direction. This is the “Achilles heel” of 
laminated advanced composites since, on occasion, loads may get rather high in this direction. 

Another problem with laminated FRP’s is that the failure modes are multiple and complex. 

These modes rely on both true material properties as well as structural properties. A homogeneous 
isotropic material fails by a single flaw, propagating to failure, whereas composites fail as a result 
of a statistical accumulation of many flaws. Unlike composites, traditional homogeneous isotropic 
materials (metals, glass, unreinforced plastics, etc.) have the luxury of enjoying a mature database 
developed over many years for predicting the effects of damage. In addition, the homogeneous nature 
of these traditional engineering materials make them much easier to work with than composites, both 
analytically and experimentally. Further complicating the problem for composites is the large number 
of ply orientations that are used in composite structures. A database must be generated for each layup 
sequence used in a structure composed of laminated composites, which can become quite cumbersome. 

Foreign object impact damage to laminated FRP is a loading condition in which all of the complex 
failure modes must be taken into account. Impact of laminated composites has been researched heavily 
in the past two decades but much remains to be understood. This research is summarized in the next 
chapter. The threat of foreign object impact to composites has hindered the use of laminates in many 
applications. Great strides have been made in developing higher strain fibers and tougher resins which 
are more damage-resistant and damage-tolerant than the early generation carbon fiber/resin systems. 
However, progress in understanding the effects of the many complex variables involved in an impact 
event needs to continue so that laminated composites can be used with even greater confidence. For the 
study presented here, continuous carbon fiber composites will be the focus, since the consequences of 
impact damage to these advanced composites pose the greatest threat. In addition, most primary load- 
bearing applications of advanced composites use continuous carbon fiber as the reinforcement in 
a polymer matrix resin, thus the emphasis on this class of material. 

In chapter 2 a review of past studies on impact damage to composite laminates is presented. 

Most of the cited literature is from experimental studies so that a general consensus of the effects of 
foreign object impact damage to composite laminates can be made. This chapter is rather broad in range 
and is intended to give the reader a good basis for understanding many of the issues involved concerning 
impact to composite laminates. Chapter 3 presents a simple analytical assessment of the effects of in-plane 
loads to change impact responses of orthotropic plates with clamped-clamped/free-ffee boundary condi- 
tions. An energy-balance solution, along with a finite element assessment, is presented. Chapter 4 contains 
the experimental procedure used in this study and chapter 5 presents the results. A discussion of these 
results is presented in chapter 6. In chapter 7, some general conclusions from this study are presented 
along with some recommendations for future work. Finally, a summary is given in chapter 8. 



2. PREVIOUS WORK 


This chapter consists of three main parts and is intended to be a fairly thorough review of the 
general problem of foreign object impact to composites. Section 2.1 is a chronological presentation 
of some of the early work that laid the groundwork for the explosive growth of studies conducted in the 
past decade. Section 2.2 covers impact damage resistance/tolerance to laminated composites. Section 2.3 
summarizes the most important features of past research conducted to date. Section 2.4 deals specifically 
with the thrust of this work and the response of impacted laminates due to a small change in boundary 
conditions. 


2.1 History and Background 

The first experiments on impact damage to composites were performed (quite naturally) on the 
first structural composite used — wood. In the early 19th century, ballistic -type impacts were studied 
by several western European naval forces in an attempt to determine the size and speed of a projectile that 
would maximize the damage to an oak ship hull. 1 As FRP came into use, the understanding of foreign 
object impact to this class of material became critical. The first experiments in the 1940’s were ballistic- 
resistance tests since the major FRP of that time, glass reinforced polyester, was being used as secondary 
load-bearing components for military aircraft.^ With the advent of “advanced composites” came the need 
to better understand the damage resistance due to a larger variety of types of impact, not just ballistic 
impacts. 

Some of the first research into the effects of foreign object impact on advanced composites 
focused on the energy needed to inflict a certain level of damage, usually fracture or penetration. In the 
mid-1960’s, glass fiber woven into cloth form, impregnated with polyester resin, was tested to determine 
the most efficient configuration to resist ballistic impacts.^ Lower velocity impact testing was originally 
performed using Charpy or Izod impact test machines, using beams made of laminated composites 
to determine what type of composite possessed the best damage resistance and what variables were 
responsible for energy absorption. 4 This was a natural extension of the type of impact testing being 
conducted on metallic materials to determine a ductility index for brittleness. With this impact technique, 
the relative importance of the fibers, matrix, and interfaces on impact resistance and energy absorption 
could be determined. In using this type of test, it was found that the fiber strain to failure was of primary 
importance in determining a composite beam’s impact resistance.^ 

As impact testing of composite laminates matured, test methods and analyses became more well 
defined. In particular, for test methods, instrumented impact and noncatastrophic testing became widely 
used and the analyses were broken down into dynamic response and quasi-static response impacts. 

2.1.1 Instrumented Impact 

By placing strain gauges, load cells, and/or accelerometers on the striker, load time, load deflection, 
and other measurements of the dynamic impact event can be made using this instrumentation. Instru- 
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mented impact testing has proven to be a valuable research tool in studying foreign object impact of 
composite laminates. The first instrumented impact tests were conducted by the British Welding Research 
Association to access the impact behavior of welded steel joints. 6 Toland was among the first to use 
instrumented impactor techniques to study the progression of damage in a composite laminate beam. 7 
By utilizing the Charpy test, the importance in distinguishing between initial elastic failure energy and 
energy absorbed during the resulting failure process was noted. About the same time (early 1970’s), 
Phillips and Tetelman published a paper demonstrating through instrumented Charpy techniques that 
linear elastic fracture mechanics could not be relied on for composite flaw analysis as it could in the 
analysis of homogeneous materials. 8 

2.1.2 Noncatastrophic Impact Testing 

Through experience gained with composite laminates, it became evident that a major threat 
is damage that cannot be seen from the outside surface. Most of this damage is due to a foreign object 
impact event. Being a laminated material, the high shear stresses set up during an impact event can cause 
the layers to delaminate in the vicinity of the impact site. In addition, with the advent of the first generation 
carbon fibers (extremely brittle) as a reinforcement, FRP laminates no longer had the impact resistance 
of glass-reinforced laminates. Thus, noncatastrophic impact testing became increasingly prevalent in the 
mid-1970’s. It is this type of testing that separated composite material impact studies from techniques 
used on metallic materials. Some of the first work involving noncatastrophic impact was performed 
on sandwich panels by Oplinger and Slepetz 9 in which a steel ball was dropped onto panels with carbon/ 
epoxy and glass/epoxy facesheets. The results clearly demonstrated the lower impact resistance 
of the carbon/epoxy facesheets. Other papers presenting results of the study of noncatastrophic impact 
to composites are found in reference. 10 Of particular interest is a study that simulated turbine blade 
impacts with steel spheres (to represent rocks), ice (to represent hail), and gelatin (presumably to simulate 
birds). 11 This study classified the failure modes produced as transverse (matrix delamination and crack- 
ing), penetrative, or breakage of the specimen. When the damage was of the transverse type (noncatas- 
trophic), the amount of damage was assessed using ultrasonic scans and cross-sectional photomicroscopy, 
two techniques used extensively today in the study of impact to composite laminates. In another study 
Husman et al. 12 examined the postimpact strength-carrying capabilities (residual strength) of impacted 
composite laminate tensile specimens. From the experimental data it was observed that the residual 
strength does drop as the impactor kinetic energy is increased and that glass fiber composites retained 
more of their initial undamaged strength than carbon fiber composites for a given impact level. Residual 
strength of impacted composites will be discussed further in section 2.3.3. 

2.1.3 Dynamic Analysis 

Since an impact event is time-dependent, that is, a load inflicted over a short period of time, 
it seems natural to use a time-dependent solution. This involves vibrations and stress waves set up 
in the material due to the impact event. This dynamic aspect of impact was examined, mostly analytically, 
fairly early among the impact studies. Moon 13 used a time-dependent solution to determine the stress 
wave propagation due to impact in anisotropic fiber composite plates. The slowest (as well as fastest) 
stress waves were associated with a unidirectional layup and had a speed of slightly less than 2 mm/psec 
(0.08 in/psec), (the fastest wave was predicted to be almost 10 mm/psec (0.4 in/psec)). Thus, if the 
impact is of a short duration, the stress waves will not yet have reflected from the target boundaries 
and will become the driving force in determining the resulting damage. Mortimer et al. 14 experimentally 


4 



measured the wave velocity in edge-impacted, cross-ply laminates as 7.38 mm/p.sec (0.3 in/psec), 
so the order of magnitude is that predicted by Moon. A previously mentioned reference 1 1 shows the 
duration of impact for a high-velocity event, on the order of 100 psec. This implies that for a 20-cm 2 
(7.9-in 2 ) specimen impacted at its center, the impacting projectile will be moving away from the target’s 
surface just as the reflected stress wave returns to the center of the plate, assuming a wave velocity of 
2 mm/psec (0.08 in/psec). In contrast, a typical low-velocity drop-weight test will be on the duration 
of 10,000 psec which gives ample time for the stress waves to decay away before the striker leaves 
the target 

Inertial effects, on the other hand, become important at much lower velocities than stress wave 
propagation effects. Inertial effects pertain to the energy needed to accelerate the out-of-plane deformations 
of the impacted plate’s mass. During low-velocity/high-mass impact, the energy to do this is negligible 
compared to the strain energies involved, but at higher velocities, the plate’s kinetic energy becomes more 
significant. Sun et al. 15 developed a continuum theory for the dynamic behavior of layered composites 
which took into account the kinetic energy of the continuum. Greszczuk 1 ^ incorporated the mass of the 
target in his analytical solution of the impact of a sphere on an elastic half-plane. Results showed that 
a higher target mass will give rise to higher impact forces and contact deformations (deformation of the 
upper surface of the target right below the impactor). Global deflections were not accounted for in this 
model. For thick laminates, it has been shown that global deflections can be ignored since they are so 
small and the material response is governed almost solely by contact stresses. 17 Inertial effects are also 
important in that the inertia of the plate and projectile can interact and cause multiple impacts. 

2.1.4 Quasi-Static Analysis 

For low- velocity impacts, many of the analyses are done as quasi-static, ignoring all time- 
dependent behavior of the event and inertial effects of the material. In essence, this analysis is that 
for a “static indentation” case in which the plate is loaded in increments of force over an essentially 
infini te period of time. Practically, this is carried out in a material loading frame. In one of the early 
studies, Oplinger and Slepetz 9 performed both low-velocity drop-weight and static indentation tests 
on graphite/epoxy and glass/epoxy sandwich panels and found that the resulting damage was similar. 
Similar results were obtained more recently by Kaczmarek and Maison 1 8 for quasi-isotropic and 
bidirectional carbon/epoxy plates. 

In the analytical approach discussed in reference 19, it was seen that a quasi-static event is defined 
by the lowest natural frequency ratio of the impactor and target, regardless of their mass ratios and impac- 
tor velocity. However, for most practical purposes, high-velocity projectiles are small objects with high 
natural frequencies which will produce a dynamic response; low-velocity impacts are really of interest 
only when the impactor mass is large and its natural frequency is low, which will produce a quasi-static 
response. 

Swanson 20 has developed a procedure based on the impactor mass and the “equivalent lumped 
mass” of the target which is not inconsistent with the work in reference 19. He showed that if the impactor 
mass is more than 10 times greater than the target’s “lumped mass,” then the impact event can be 
considered quasi-static. For various plates and beams of different boundary conditions, the “lumped 
mass” was seen to be no more than one-half the entire plate or beam mass. Thus a crude, but effective, 
rule of thumb is that the impactor must be more than 20 times the target mass for the impact event to be 
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considered quasi-static. However, boundary conditions can affect this rule if the specimen is supported 
such that it is stiff in flexure. Highsmith 21 found that quasi-static results produced slightly more damage 
than impact events of similar maximum load. The specimens in his study were 20 plies thick and clamped 
over a 6.35-cm- (2.5-in-) diameter opening. 

Jackson and Poe 22 demonstrated that impact force alone can be used to quantify an impact event 
for a given material and layup but only when the impactor mass is large and the velocity is low. In their 
study it was found that the maximum delamination damage area was a linear function of the force 
of impact whether the force was introduced by falling weight or static indentation tests. Furthermore, 
the support diameter of the impacted plate had no effect on the delamination size when impact force 
was used to quantify the event. However, for large plates, deviation from quasi-static behavior was noted 
since the stress waves set up by the impact event took longer to travel the length of the specimen and 
reflect back from the boundaries, resulting in some dynamic effects. 

2.2 Damage Resistance and Damage Tolerance 

Damage resistance and damage tolerance are two distinct and completely different aspects 
of impact damage to laminated composites. Damage resistance refers to the composite material’s ability 
to minimize matrix cracking, delamination, and fiber breakage due to an impact event. Damage tolerance 
refers to a material’s ability to retain its desired properties once damage has been incurred. This is usually 
referred to as “residual” properties. 

During an impact event, many parameters determine a material’s damage resistance. Structural 
parameters such as boundary conditions and the stress state during impact can greatly affect a laminated 
composite’s damage resistance. It is the intent of this study to assess the damage resistance effects of 
a tensile stress state being applied to the target during an impact event. This will usually be referred to as 
a “tensile preload” throughout the remainder of this paper. Tensile prestress and tensile prestrain will also 
be used but implies the same physical phenomena. 

2.3 Key Findings of Past Research 

The following review concerns only fiber-reinforced laminated plastics being transversely 
impacted by a foreign object. The types of damage incurred during an impact event will be examined 
and the effect that the stacking sequence has on this damage will be investigated. The next extension 
of the impact problem involves residual load-carrying capabilities of the damaged laminate. A discussion 
of how different fiber/resin systems and boundary conditions can affect the amount and type of damage 
induced in an impacted composite will then be examined. Finally, a very important topic to those working 
with actual structures made of reinforced plastics, seating from laboratory specimens to larger entities, 
will be discussed. 


6 



2.3.1 Damage Modes 


There are three major modes of damage that can be incurred during foreign object impact events. 
These are matrix cracking, including fiber/matrix debonding; delamination; and fiber breakage. The 
formation and extent of these damage modes are heavily dependent upon specimen geometry, or more 
precisely, the ability of the impacted structure to deform. For relatively stiff structures, contact forces will 
dominate, and for flexible structures, in-plane forces will dominate due to global bending. The contact 
force generated by the impactor, if sufficiently large, can cause transverse matrix failure in the impacted 
ply and/or those beneath it parallel to the direction of fibers in each ply, due to the high shear stresses, x yz , 
that are set up (shown schematically in figure 1). 



Figure 1. Top ply cracking due to contact stresses. 

As the impactor drives itself into the specimen, the material from the top ply tends to spread 
the load over a wider diameter, producing matrix cracks slightly away from the impact site in die next-to- 
the-top ply. As these transverse cracks develop, they become wider apart through the thickness as shown 
schematically in figure 2. These transverse cracks usually join one another through the acute angle with the 
adjacent ply via delamination, shown schematically in figure 3. When cross-sectioned, the transverse 
cracks and delamination will form a characteristic cone shape as noted in references 23-27 (shown 
schematically in figure 4). The through-the-thickness assimilation of the delaminations will give the 
characteristic “stair-step” shape, noted in references 28 and 29. 

For composite laminates that experience large deformations during an impact event, large in-plane 
stresses can cause damage to form. The intrinsic weakness of a ply perpendicular to its fibers causes 
the bottom plies of heavily deformed laminates to split along the fibers in that ply, due to the high tensile 
strains from bending. The transverse shear in the plate will be at maximum near the center of the plate 
and it has been observed that this region contains the planes of largest delaminations, except when spalling 
occurs on the nonimpacted face, after an impact event. 29 
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+ 45 ° 


0 ° (Splits parallel to direction of cut and thus cannot be seen in 0° plies) 



0 ° 


Figure 4. Schematic representation of cross-sectional view of transverse cracks 

and delaminations due to contact forces. 


By using simple bidirectional layups (fibers aligned only in the 0° and 90° directions) and 
examining the resulting matrix cracking and delamination, more insight is given into the problem 
of damage development in impacted composite laminates. In tests by Chow et al ., 30 for [0„/90 16 _2 n /0 n ] 
(n=6 and 7 ) stacking sequences, damage was seen to form as matrix cracks in the 90° plies with extensive 
delamination between the 0° and 90° plies (when sectioned parallel to the 0° plies). Smaller matrix cracks 
were observed in the 90° plies that did not join with the delamination, thus all matrix cracks do not end in 
delamination propagation. In these tests, delamination was observed directly under the impact site, 
in dicatin g that flexural behavior is responsible for this growth. These experiments were unique in that the 
laminates were supported as wide beams (plates clamped-clamped/free-free) and impacted by a line nose 
impactor rather than the conventional spherical-shaped nose impactor. Thus, cracking in the 0° plies was 
limited since bending and deformation were only longitudinal to the fibers in these plies. A spherical- 
shaped nose impactor typically causes the laminate to deform in a “bowl” shape, causing matrix splitting 
between fibers in all of the plies. By dispersing the plies more in a [O 3 , 9 O 3 , 03 , 9 O 3 , 03 , 9 O 3 , O 3 ] layup, it 
was found that matrix cracking still occurred in all of the 90° plies but the delaminations between the 0° 
and 90° plies became progressively smaller from bottom to top. 

A finite element analysis of delamination growth from matrix cracks was given in reference 31. 

It was demonstrated that normal stresses develop at the matrix cracks that are tensile on the upper interface 
toward the impacted side and tensile on the lower face away from the impact site. Thus, delaminations are 
predicted to grow toward the impact site on the upper interface and away from the impact site on the lower 
face, producing larger delaminations near the bottom plies when bending is the governing stress producer. 

These findings were verified for quasi-static loading of bidirectional laminates in three-point bends 
while observing a polished edge in reference 32. 

An experimental microstructural analysis of impacted [ 0 / 90 / 04 ] s and [0/90/0 revealed that 
the more 0/90 interfaces in a laminate, the smaller the delaminations . 33 In addition, delamination growth 
appeared to grow only away from the impact site on the bottom interface of the 0/90 boundary. No trans- 



verse cracking was observed in the 90° plies and the delaminations were observed only at every other 
layer. The authors attributed this to the processing of the panel, but as demonstrated in figure 5, this 
is due to the nature of how the specimen was sectioned coupled with how the delaminations form 
as shown in figures 1-3. 

For high-velocity impact loading, arguments have been given that the tensile stress waves 
are responsible for matrix cracking on the upper surface of an impacted composite. 34 Evidence 
was also given that matrix cracking and delaminations interact simultaneously since matrix cracking 
was not observed in delaminated areas in interfaces away from the impacted side of a bidirectional 
composite. 
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Figure 5. Schematic of cross-sectional view of bidirectional laminates. 


The process of delamination formation being the result of bending stiffness mismatch between 
plies of different orientation was given by Liu. 35 He proposed a bending stiffness mismatching coefficient 
that would qualitatively give the size and shape of the delamination, resulting from the differences in the 
L >1 j terms of adjacent plies. The shape was seen to be of a “peanut” type, mathematically a lemniscate. 
The size was seen to increase as the angle between plies approached 90°. However, a host of empirical 
data 28 ’ 29 * 36-42 shows that the characteristic shape of delaminations between plies is mostly governed 
by matrix cracks and the lemn iscate has a significant magnitude in the 45° direction, whereas it is 
predicted by Liu to be zero in this direction. 

Fiber failure occurs in laminates experience seeing such large contact stresses that surface fibers 
are broken or such large bending stresses as to exceed the strength of a lamina in the direction of the 
fibers. A detailed discussion with empirical data for fiber breakage in a laminate due to large deflections 
is given in reference 43. This study showed that the onset of fiber breakage was in laminae near 
the bottom of the plate. 


10 



2.3.2 Stacking Sequence 

Stacking sequence of the laminate is very important in determining the mode and size of damage 
due to foreign object impact; thus, the residual performance of the structure. For structures which 
experience high contact stresses, it has been shown 44 that placing nonload-bearing fibers on the outer plies 
of the laminate to protect the load-bearing 0° plies is advisable. However, as Cantwell et al. 45 suggested 
outer 0° plies would increase the elastic energy absorption capacity by increasing the flexural stiffness. 
When examined from a damage-tolerance point of view, Cantwell et al. 45 showed that the residual com- 
pressive strength of a composite with the load-bearing 0° fibers away from the surfaces and not next 
to the outermost delamination was larger than specimens that had outer 0° fibers that delaminated 
from the angled ply next to them. The delaminated 0° outer plies could easily buckle away from the rest 
of the laminate and the strength effect of these fibers would then be lost. 

From a damage-resistance point of view, Kwon and Sankar 27 examined 0/90 bidirectional and ji/4 
and jt/8 quasi-isotropic stacking sequences. Results showed that initial damage due to transverse loading 
occurred relatively early in the jt/8 laminates but this damage did not cause as drastic a load drop 
(as measured by instrumented impact) as the other two layups. This is because transverse cracks between 
plies of different orientation have a smaller angle between them in the Jt/8 laminates and the resulting 
delaminations that form between these cracks between plies (see section 2.3.1) sweeps out a smaller area 
which requires less energy. The bidirectional laminates have a full 90° between ply orientations and the 
delaminations between cracks in adjacent plies in these laminates must sweep out the largest area of any 
stacking sequence, therefore, the large load drop due to the formation of these large delaminations. In this 
study, 27 cross-sectional examination of the damaged laminates revealed that the bidirectional laminates 
contained the most massive and numerous delaminations followed by the ji/4 and then the Jt/8 laminates. 

By stacking a laminate such that the number of grouped plies of equal orientation is at a minimum, 
increases the interfaces at which a delamination can occur (delaminations do not form between plies 
of equal orientation). This was shown to reduce the overall damage area in references 46 and 38 since 
the laminates with less ply interfaces at which delamination could occur produced larger delaminations 
between the plies that could delaminate in order to dissipate the same amount of energy. For bidirectional 
laminates it has been shown 47 that dispersing the plies of different orientation significantly increases 
a laminate’s impact resistance (as measured by impact energy to produce initial damage) and that 
increasing the number of 0° plies on the outer surfaces increases a laminate’s compressive damage 
tolerance by increasing the buckling resistance of the sublaminates on the critical outer plies. 4 ** 

For honeycomb sandwich panels, 4 ^ it was found that the facings would experience less delamina- 
tion than equivalent laminates and that contact forces dominate since bending is restricted on the facings 
of sandwich panels due to the existence of the core. This was supported by data which showed that at low 
impact, energies that did not cause large laminate deflections (and were thus contact stress-dominated), 
facings and equivalent laminates had approximately the same total delamination areas, but as the impact 
level increased to cause large deflections of the laminates, the facings would experience a smaller, total 
delamination area since the large bending stresses in the laminate caused larger delaminations than those 
caused by contact stresses. It was also found in this study that layups with large angular orientations 
between plies (e.g., 90°) tended to cause a larger total delamination area across all of the ply interfaces 
(through-thickness) compared to layups that had 45° as the largest angular orientation between plies. 
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2.3.3 Residual Tensile Strength 


In this and the following section, the effects of incurred impact damage on composite laminates 
is assessed by examining residual strength (tension or compression), a measurement of the composite’s 
damage tolerance. While some papers will present results as plots of residual strength versus damage area 
(thus isolating the damage tolerance response), many others show results on residual strength versus 
impact energy plots which actually serves to effectively combine the impact resistance and damage 
tolerance portions of the impact event. 

Once damage has been introduced into a composite laminate, the critical question is how the dam- 
age will affect the desired properties of the material. Measurements and predictions of residual strength 
have seen extensive research, especially since the early 1980’s. Much of this work has concentrated on 
compression-after-impact (CAI) strength since impact events cause delaminations which can lead to a 
lower compressive strength due to the localized instability caused by the delaminations (this will be dis- 
cussed in the next section). However, many load-bearing composite structures will be dominated by ten- 
sile forces (such as pressure vessels), so this aspect of the impact problem also needs to be examined. 

Husman et al. 12 were among the first to investigate and attempt prediction of the residual tensile 
strength of impacted coupons. The residual strength prediction was based on an analogy with an inserted 
flaw of known dimensions and a stress intensity factor for an isotropic material. As a result it was pre- 
dicted that the residual strength would drop off in proportion to the square root of the impact energy 
of the projectile. It was noted that departure from this analysis will occur as the specimen experiences 
complete penetration. Graphically, the residual tensile strength versus impact energy curve as predicted 
by Husman et al. is given in figure 6. Limited data were given to support this analysis and since ballistic- 
type impacts were used, the fibers on the impacted surface (which were 0° plies) may have broken at low 
impact energies (but still high velocities), yielding a drop in residual at the lowest impact energies used 
in this investigation. 
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Figure 6. Residual strength versus impact energy curve as predicted by Husman et al. 12 

This study also suggested that the width of the specimen be at least six times the diameter of the 
projectile in order to minimize damage zone/edge interaction. This is an important factor to consider when 
evaluating residual strength on laboratory specimens. For example, one researcher may impact 2.54-cm- 
wide (1-in- wide) specimens while another performs identical impact tests on 15.24-cm- (6-in-) wide 
specimens. If the damage zone is 2.54 cm (1 in) in both specimens, then the narrower specimens will 
register a much lower residual strength (possibly zero; i.e., catastrophic failure) than the wider specimens 
which possess undamaged material around the damage zone to help carry the load. Thus, if a structure 
is larger than the coupons used to perform the impact damage characterization on it, then the coupons need 
to be of sufficient width so that free-edge stress/damage zone interactions do not occur. In addition, 
if the damage zone is a sizable percentage of the specimen width, then the cross-sectional area rendered 
ineffective due to damage must be removed from the residual stress calculations. This is an argument for 
measuring residual strain which eliminates needing to know the cross-sectional area rendered ineffective 
to carry load. 

Data by Rhodes, 50 also generated by ballistic-type impacts, seems to support the notion that both 
residual tensile and compressive strengths drop rapidly from the undamaged strength even at the lowest 
impact energies tested (fig. 7). It is worth noting that these data were from honeycomb sandwich panels 
whose composite facesheets undergo smaller deflections (due to the backing rigidity of the core), resulting 
in a more localized damage which includes more fiber breakage. 
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Impact Energy 


Figure 7. Rapid reduction in residual strength as measured by Rhodes. 50 

Other studies have shown that for high-velocity impacts, the residual tensile strength does indeed 
begin to drop quite rapidly at the lowest impact energy levels tested. 51-54 A study by Cantwell and 
Morton 26 contained residual tensile strength data for eight-ply specimens in which the tensile strength did 
not begin to drop off until a critical impactor kinetic energy was obtained. Cairns 55 also found that there 
existed a plateau region of residual strength before degradation for 12-ply carbon/epoxy specimens. This 
study also showed that for equivalent energy, higher velocity (smaller mass) projectiles, the residual 
strength began dropping off much more rapidly, indicating that the plate’s inertial effects become more 
important at higher velocities, as pointed out in section 2.1.3. 

Lower velocity impacts, such as drop-weight tests, generate a different material response than 
high-velocity impacts, as mentioned in sections 2.1.3 and 2.1.4. Therefore, it may be expected that the 
residual tensile strength response may differ. Comparisons of equal energy high- and low- velocity 
impacts have been conducted by Adsit and Waszczak 56 and Dorey. 44 * 57 In reference 56, carbon/epoxy 
facesheet honeycomb panels were impacted at low velocities with drop weights and at high velocities 
utilizing a compressed air gun to project small granite stones. A plot of residual tensile strength versus 
impactor kinetic energy shows that the strength reduction begins at a lower level of impact energy for the 
high-velocity stones and also drops off at a faster rate with increasing impactor kinetic energy for the 
higher velocity tests. In reference 44, the residual strength is seen to drop at a smaller impact energy for 
ballgun impacted specimens, but once a decrease in strength is seen, it is more rapid for drop- weight tests. 
This trend is shown schematically in figure 8. This effect is much more pronounced for residual flexural 
testing as seen in reference 56. The type of resin used can have an effect on these results, as discussed 
in section 2.3.5. 
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Figure 8. Schematic of residual tensile strength versus impact energy 
for high- and low- velocity impacts. 

Low-velocity impact tests with resulting residual tensile strength tests show that, for most cases, 
a region of impact energies where strength degradation does not occur is readily apparent. This 
phenomena is seen in data by Cantwell et al., 45,5 ® porey et al., 5 ^ Morton and Godwin, 60 Caprino,^ 
Nettles, 62 Bishop, 63 Lai, 64 and Poe. 65 

However, once strength degradation does occur, the residual strength usually decreases rapidly 
with increasing impactor kinetic energy, depending upon the material and geometry of the impact event. 
For example, specimens supported as short beams 60 " 61 or as small diaphragms 62 experience a high level 
of contact stresses immediately below the impactor, thus breaking fibers and causing a rapid reduction 
in residual tensile strength. Specimens supported over larger areas, such as in references 59 and 64, 
demonstrate a more gradual decrease in residual tensile strength with increasing kinetic energy since 
the large elastic deformations can absorb much of the incident impactor kinetic energy and local contact 
forces under the impactor are less severe. The effects of boundary conditions will be examined further 
in section 2.3.6. 

2.3.4 Residual Compressive Strength 

Residual compression strength after impact has drawn much more interest than residual tensile 
strength since laminates have been shown to be much more susceptible to this type of loading. The impact 
event can cause delaminations within the material which cannot be seen from the surface, yet can severely 
reduce the CAI strength of the composite structure. While the tensile strength of a composite laminate is 
fairly dependent on the test sample’s geometry, compression testing of laminates is virtually a complete 
structural test as opposed to a true material test. Thus extreme caution must be used when examining data 
from CAI tests and direct comparisons between studies is virtually impossible. Therefore, only general 
observations will be made based on the results of some of the findings of past research on this topic. 
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In low-velocity impact studies where both residual tensile and compressive strengths were 
tested, the residual compressive strength always fell at a faster rate than the residual tensile 
strength. 44 * 43 * 57, 60, 63, 66 Rhodes^ 0 presented data for projectile-type impacts that showed the residual 
tensile strength decreasing at a slightly faster rate than the residual compressive strength. These results 
are because a specimen will lose tensile strength only when load-bearing fibers are broken, but com- 
pressive strength can drop due to delaminations. In most test programs the specimens are supported 
such that matrix cracking and delamination will occur before any fibers are broken. Broken fibers without 
delaminations only occur in thick specimens, or in ballistic-type impacts, in which the specimen does not 
have time to deform and create the high interlaminar shear stresses that cause delamination. 

For low- velocity impacts (drop weight), the characteristic shape of the CAI versus impact energy 
curve is heavily dependent upon the specimen material and geometry. In a study by Prichard and Hogg, 67 
carbon/epoxy and carbon/polyetheretherketone (PEEK) specimens were supported over a relatively small, 
clamped annular opening to induce delamination-type damage. The resulting CAI versus impact energy 
curves show that the compression strength of the carbon/epoxy dropped off extremely rapidly whereas the 
curve for carbon/PEEK gradually dropped off with increasing impactor kinetic energy, almost in a linear 
fashion. The same results were found in references 60 and 68 for specimens clamped with a larger annular 
diameter support, although the carbon/epoxy CAI strength did not drop off quite as rapidly. In general, it 
has been shown that for a wide range of carbon/epoxy specimen geometry and CAI test methods, the CAI 
strength falls rapidly for increasing impactor kinetic energies at the lower level of impact energies until 
a plateau or modest linear decrease zone is reached. 44 * 45 * 61 , 63 , 66 , 69-75 


For high-velocity- (projectile-) type impacts the CAI strength versus impact kinetic energy curve 
does not possess as drastic a drop in CAI strength with increasing impactor kinetic energy as shown 
in references 50, 52, 76, and 77. This is mainly because the much more localized response of the high- 
velocity impacted specimens in which delamination tends to be extensive through the thickness, but not 
as widespread in areas such as lower velocity-impacted specimens. Thus, the lower velocity-impacted 
specimens contain a larger delamination area, causing a larger reduction in structural stability of the test 
coupon than the smaller delaminations formed by higher velocity impacts. This notion is supported 
by data in reference 78 in which the CAI strength of specimens impacted with high- and low- velocity 
impactors did not vary a great deal when plotted versus damage width, but were different when plotted 
versus impact energy. 

From drop-weight impact data, Prichard and Hogg 67 showed that when plotted as a function 
of damage width, the residual compression strength of two very different laminated composite materials 
(epoxy and PEEK matrices) were the same. Similar results were obtained by Teh and Morton 68 and 
Dempsey and Horton. 73 However, other studies 79 have shown that a specimen of one type of material 
may have a much larger damage zone than another type of material impacted under identical conditions, 
yet have a higher retention of initial compressive strength. Indeed, if the residual strength versus damage 
size plots were similar for all materials, then this would imply that there is no difference in the damage 
tolerance of laminated composites, only a difference in damage resistance. 

Of practical importance for composite structures that see multiple impacts at nearby sites (such 
as during a hailstorm), Jones 80 showed that the CAI strength was unaffected by the multiple impacts 
as long as they were approximately one damage zone diameter apart (i.e., the impact sites do not interact 
when spaced apart by at least one diameter of damage zone). 
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2.3.5 Fiber/Resin Systems 

Of the most interest to designers and end users of composite structures are which fiber/resin 
systems are the most damage resistant, or more importantly, damage tolerant. Since the chosen fiber/ 
resin system is a variable that can be controlled from the onset of design, this information is vital to any 
program utilizing laminated composites. In the past decade, a concentrated effort has taken place within 
the industry to produce fibers with higher strain to failures and matrix resins that better resist cracking and 
delaminations. More exotic techniques such as stitching or using three-dimensional preforms, 81-85 
interleaving, 86-88 and using hybrids (more than one fiber type) 36 ’ 89-95 have been investigated, although 
they have not yet seen widespread usage. 

The toughness of the matrix resin can severely affect the damage resistance/tolerance of laminated 
structures. A direct comparison of thermoplastic and thermoset resins usually shows that the thermoplas- 
tic matrix material is more damage resistant. The most promising high-performance thermoplastic resin 
is PEEK which has been studied extensively in damage-tolerance programs. A comparison of PEEK 
with epoxy resin composites was made by Dorey et al. 59 and Demuts and Sharpe 96 and showed that the 
PEEK laminate was much more damage resistant (as measured by C-scanned damage area), resulting in 
greater retention of CAI strength. In reference 59, for the highest impact energies used (8 J), the carbon 
PEEK had about one-fourth of the damage area as the carbon/epoxy tested. When compared to another 
thermoplastic such as polyphenylene sulfide (PPS), 97 PEEK was still more damage resistant with about 
a 50-percent smaller delamination size for the highest impact energies tested. 

By plotting CAI strength versus damage size for carbon/epoxy and carbon/PEEK laminates 
and obtaining similar results, 67 it follows that the superior CAI resistance of the PEEK laminates is due 
mainly to the smaller delamination area formed for a given impact energy level. This is further supported 
by work performed on a number of toughened epoxies and bismaleimide (BMI) resins in which the CAI 
strength was found to correlate strongly with mode II interlaminar fracture toughness. 98 Mode II fracture 
resistance controls the amount of delamination that forms during an impact event whereas the CAI 
strength is strongly controlled by mode I peeling as the delamination grows. 

Although thermoplastics such as PEEK have distinct advantages over thermoset resins, these 
materials require a high temperature and pressure to cure the plies into a laminate. Furthermore, lack 
of any tack and little drape make further processing complicated. Until more advanced processing 
methods are found, thermosets will continue to dominate the resins in use on composite structures; thus, 
mechanisms for making existing thermosets tougher has seen much activity. The two principle methods 
of toughening thermosets is by modifying these resins either with elastomeric particulates or by blending 
in a controlled amount of thermoplastic resin. Both methods can substantially increase a laminate’s 
toughness and modified resins of both type are commercially available and in wide use. Damage 
resistance/tolerance comparisons of some of these resins have been made in a number of 
studies. 53, 75, 79, 99 
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2.3.6 Boundary Conditions 

Another critical factor that determines the amount of damage that a composite laminate will 
experience due to a foreign object impact event is the laminate’s geometry and boundary conditions. 

For high-velocity impacts, the boundary conditions of the target are less important than the areal mass 
of the target in determining the impact response 100 (unless stress waves account for damage). The 
impactor has already rebounded from, or penetrated the target by the time the stress response reaches 
the edges of the target specimen; thus, boundary conditions are much more critical for lower velocity 
impacts. Cantwell and Morton^ 0 ' showed that for a given target thickness, the increase in energy required 
to cause initial damage and to cause perforation for impacted beams, was roughly proportional to the 
increase in beam length. It was also shown in this study, that for a given impact energy, lengthening 
the beam caused the beam to experience larger deflections, thus storing more of the energy as elastic 
energy and reducing the size of the damage zones. Kwon and Sankar 27 demonstrated that if the 
delamination size is plotted as a function of maximum load of impact, then the specimen support size 
effect is eliminated, since larger supports allow larger deflections and thus lower impact forces. In a study 
by Verpoest et al ., 102 the size of the support was shown to be more important than the shape. Another 
aspect of the varying response of the target, due to support geometry, is where the majority of damage 
is most likely to form. For specimens supported in a smaller geometry, the top surface will experience 
higher contact stresses and thus more damage will form at this location than on a similar target, only 
with larger geometrical boundaries, which will experience bottom surface splitting/fiber breakage due 
to the larger bending stresses that are present . 103 

Also of importance is how the specimen is supported at its boundaries (i.e., clamped, simply 
supported, etc.). The same general rule of thumb applies here also. If the boundary conditions restrict 
the specimen from deforming, then the damaging stresses will migrate more from bending type to contact 
type. This was shown analytically by Shivakumar et al. 10 ^ 

Another critical aspect of the impact event, the shape of the impactor, has seen relatively little 
attention. Various sizes of hemispherically ended impactors have been used, but other geometries 
of impactor were studied by Poe . 105 This study was conducted on thick sections of graphite epoxy that 
simulated a section of a rocket motor case. Impactors of two different hemispherically ended radii — a 
sharp, bolt-like rod — were used. The comer-type impactor caused very high localized Hertzian contact 
pressures, producing visible damage earlier than the hemispherical indentors. The cylindrical bolt-type 
indentor acted like a punch, slicing through the upper layers of the laminate. The shape of the impactor 
will have a much greater effect on laminates which are relatively stiff, since the Hertzian contact pressures 
are more important for this type of impact event. 

2.3.7 Scaling of Impact Damage 

The ultimate goal of the vast amount of research that has been conducted on impact-to-laminated 
advanced composites is to be able to predict the effects of a foreign object impact event to a full-scale 
composite structure. Since this can turn into a very expensive program, laboratory-size specimens must 
be used and the results scaled to those of larger structures. The problem of scaling impact damage 
to composites has been investigated by Morton . 106 Morton used the principles of similitude and pointed 
out, that for most practical cases, complete similarity between “model” (coupon) and “prototype” (full 
scale) is rarely possible when rate- or notch-sensitive materials are involved. For composite laminates, 
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the problem is further complicated since the fiber diameter and lamina. thickness should ideally be variable 
to accommodate scaling. In practical applications, this is virtually impossible. Test results on impact 
of composite beams showed that classical scaling laws apply for elastic behavior of the beams (i.e., before 
damage is induced) but significant deviation was observed as damage was induced into the specimen with 
the smaller specimens always stronger than the scaled-up ones. Swanson 1 has also examined the scaling 
effects of impact damage to composite laminates and found that delamination deunage was dependent 
upon the absolute specimen size, whereas, fiber breakage was not. Thus, the effects of seeding must be 
different for different failure modes. 

For delamination damage, Bucinell et al. 10 ^ tested laminated plates of scaling factors 1, 2, and 4 
emd showed that the response followed an inverse square root of the scaling factor. The response mea- 
sured was the percentage of laminate area delaminated as detected by ultrasonic C-scans. The impact 
parameters and responses were scaled according to the development of Morton. 106 It was also demon- 
strated in this study that quasi-static and dynamic impact events of equal energy scale similarly. 

Through finite element analysis, Chen et al. 109 showed that the stress concentration due to 
a damage zone is concentrated around that zone and local stress distributions between test coupons and 
stringer-stiffened panels were similar. Thus, the measured differences in CAI strength in which the larger 
panels retained more strength 1 10 is due to the larger panels experiencing less damage for a given impact 
energy per unit thickness. 


2.4 The Effects of In-Plane Loading, N x 

A limited amount of research into the effects of tensile or compressive preloads on the impact 
response of composites has been performed. Olster and Roy 111 performed ballistic penetration tests on 
tensile preloaded graphite/epoxy plates and found that the combination of impact and preload could cause 
a 45-percent reduction in ultimate tensile strength whereas penetrated plates without a preload would fail 
at an approximately 38-percent reduction in ultimate tensile strength. This suggests an interaction between 
the preload stresses and the stresses set up by the impact event. Since all specimens were completely 
penetrated by the projectiles, the residual tensile strength was found to be independent of the preload. 

For nonpenetrating impacts on preloaded specimens, Awa (formerly Sharma) 1 ^ -1 16 has 
produced the most experimental data on this topic. Typically, a “failure threshold curve” is obtained 
by plotting specimen preload (tension or compression) versus impactor kinetic energy and drawing 
a faired curve between the data in which the specimens broke upon impact (a “catastrophic” break) 
and those which survive the impact. A sample curve (not from actual data) is presented in figure 9 
for clarity on this subject. As will subsequently be seen, a well-defined area between the two zones 
exists in most actual experimental data. 
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0 = Failed On Impact 


Q = Did Not Fail On Impact 



Kinetic Energy of Impactors 

Figure 9. Example plot of preload versus impact energy curve 
to demonstrate concept of a failure threshold curve. 

In a study on honeycomb sandwich panels impacted by projectiles from a light gas gun 1 ^ 
it was found that tensile preloads on the facings of sandwich panels could cause catastrophic failure 
at stress ratios (ratio of preload stress to ultimate breaking stress) below 40 percent for 1 J (0.7 ft-lb) 
of impactor kinetic energy. This was for facings of T300/5208 carbon/epoxy in a stacking sequence 
of [90,+45,-45,0] s bonded to a 130 kg/m 3 (8.1 lbm/ft 3 ) aluminum honeycomb core. For stacking 
sequences of [+45,^45, 0 4 ] s and [+45,-45,90,0] s , higher prestresses needed to be applied to cause 
catastrophic failure for a given impact energy. This is understandable for the [+45,-45,0 4 ] s layup since 
it is thicker and contains more 0° plies. Limited data were supplied for the [+45,^45,0 4 ] s and 
[+45,-45,90,0] s specimens, so conclusions on why this may be are difficult to draw. Contact stresses 
surely dominated in this study since the facesheets were backed by a relatively hefty honeycomb 
of 130 kg/m 3 (8.1 lbm/ft 3 ) density. For specimens that did not catastrophically break upon impact, 
the residual strength was obtained. Although this aspect of the study was not highlighted in the paper, 
it is observed that the tensile residual strength is not much higher than the failure threshold curve 
for the limited data presented. For compression, the residual strength could be much more above 
the failure threshold curve. This is probably due to the stiff honeycomb backing preventing large 
deflections; thus, the large delamination areas associated with them. 

In a study that contained much more data, U 3 the effects of a preload on the residual strength 
of specimens were noted. The specimens were made of T300/5208 carbon/epoxy with a stacking 
sequence of [452,—452,02,902] s . The support conditions were clamped/clamped (perpendicular 
to the 0°-direction)-free/free with a width of 7.62 cm (3 in) and length (0°-direction) of 15.24 cm 
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(6 in) for tension loaded specimens and 6.35 cm (2.5 in) for compression-loaded specimens. An 
aluminum sphere of 1.27 cm (0.5 in) diameter was propelled at velocities between 18.6 m/sec (61 ft/sec) 
and 106.7 m/sec (350 ft/sec) and struck the specimens at their geometrical midpoints. This testing 
geometry will induce relatively large, local contact stresses while not causing much global deflection 
as would be present in a high-mass/low-velocity impact event, thus minimizing the increased stiffness 
response due to the tensile preloads (or decreased stiffness response due to compressive preloads). 

For the data presented on tensile loading, a faired curve could easily be drawn between regions 
of catastrophic and noncatastrophic failures. For the specimens that did not experience catastrophic 
failures, the residual tensile strengths were much above the catastrophic failure curve. Specimens which 
were impacted while under no preload demonstrated a residual strength versus impact energy curve higher 
than the catastrophic failure curve. The residual strength data from the preloaded specimens that did not 
catastrophically fail fall close to the no preload residual strength curve. An additional set of tests were 
conducted using the catastrophic failure curve as a guide for setting preloads at three different impact 
energies. All of these specimens were preloaded in the region of the catastrophic failure curve. As 
expected, most of the specimens failed catastrophically upon impact, and those that did not had very 
little residual strength remaining. 

Similar experiments were conducted for specimens under a compressive preload. For the 
specimens that did not fail catastrophically upon impact, the residual compression strength was seen 
not to rise much above the catastrophic failure curve. In addition, the catastrophic failure curve was close 
to the residual strength versus impact energy curve of specimens impacted with no preload whereas for 
the tensile specimens, the catastrophic curve was approximately 30-percent lower than the residual 
strength versus impact energy curve for unloaded specimens. These results are sketched out graphically in 
figure 10 for clarity and comparison. Only pertinent data points are included to reduce clutter in the plots. 
Thus, from this study it appears that a tensile preload can reduce the residual strength-carrying capability 
of the impacted laminate, but only when preloaded near the “catastrophic failure” curve. 
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Tensions Compression 

0 — Catastrophic Break — ^ 

O — Survived Impact — ^ 

Q — Residual Strength — ^ 

□ _ Residual Strength — ^ 

(No Preload) 



Kinetic Energy/Thickness (J/cm) 

Figure 10. Representative data from Avva. 113 

In a very similar study, 114 specimens of T300/934 carbon/epoxy with a stacking sequence 
of [+45,-45,0,90]2 S were impacted by an aluminum sphere of 1.27 cm (0.5 in) diameter while under 
a tensile preload. The catastrophic failure curve for this study is very similar to that shown in figure 10. 
The only noticeable difference between these two studies is that the residual strength of the preloaded 
specimens that did not fail catastrophically fell closer to the catastrophic failure curve in this study as 
compared to the study in reference 1 13. A summary of this study can be found in reference 1 15. 

In an impact damage study by Starnes et al. 76 the residual compressive strengths of 48-ply laminates 
loaded with a compressive preload were evaluated. This represents an extremely thick test coupon 
and global deflections should be negligible due to the laminates high bending stiffness. The projectile 
used was a 1.27 cm (0.5 in) aluminum sphere propelled between velocities of 52 to 101 m/sec 
(171 to 331 ft/sec). The material used was T300/5208 carbon/epoxy with a stacking sequence 
of [+45,-45,02,-1-45,-45,02, +45,^15,0,90]2s cut into specimens 1 1.4 cm (4.5 in) wide by 24.8 cm 
(9.75 in) long. Catastrophic failure threshold curves were determined from the data and the two zones 
between specimens that survived the impact and those that failed on impact were very distinguishable. 
For those specimens that did survive impact the residual strengths were determined. Below 5 J (3.7 ft-lb) 
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of impact energy, the residual strength values were within the scatter band of specimens that had not been 
impact damaged. For impact energy values that did cause a reduction in residual compression strength, 
the preload appears to decrease the residual compressive strength to a small degree, although only six data 
points are used for this conclusion. In a continuation of this study, 1 ^ additional data on 40-ply laminates 
of stacking sequence [+ 45 ,- 45 , 04 ,- 45 ,-t 45 ,+ 45 ,- 45 , 04 ,- 45 , 445 ,+ 45 ,- 45 , 02 ] s demonstrated that 
specimens impacted while under a compressive preload experienced slightly larger areas of damage, 
as detected by brittle lacquer spallation on the nonimpacted side of the test coupon. Specimens of a quasi- 
isotropic layup, [+45,-45,0,90,-45,+45,0,90 ] 3s were also tested, with and without compressive preloads, 
with no real discernible drop in CAI strength due to the compressive preloads. These results are also 
presented in reference 117. 

Low-velocity drop- weight tests on tensile prestressed thinner specimens were performed 
in a study by Park. 118 T300/914C carbon/epoxy with four- and eight-ply bidirectional layups were 
impacted by a 1.5-kg (3.3-lbm) mass impactor with a 6-mm (0.24-in) diameter tup. The 2.5-cm- (1-in-) 
wide specimens were supported as clamped-clamped/free-free with a span of 5 cm (2 in). Catastrophic 
threshold curves were generated (termed “impact maps” in this study) and showed that ballistic impacts 
had a higher threshold curve than drop- weight impacts since the ballistic tests caused localized penetration 
and the drop-weight tests caused damage across the entire specimen area. A limited number of tests on 
the residual tensile strength of preloaded coupons, impacted ballistically, showed that the residual strength 
dropped sharply as the preload neared the ultimate tensile strength of the test coupon. 

A C-scan analysis of impacted beams under a tensile preload was conducted by Sankar and 
Sun. 119 Although only one initial preload value was used (about one-third of ultimate load) the results did 
show differences in those beams impacted with an initial tensile stress and those impacted with no initial 
stress. The beams were 20-ply graphite/epoxy with a stacking sequence of [02,902, 02 , 902, 0 2 ] s supported 
as clamped-clamped/free-free. A 1.27 -cm- (0.5-in-) diameter steel ball was used as a projectile at speeds 
of 10-40 m/sec (33-131 ft/sec). The span of the impacted beams was set at 18 cm (7 in) and were 3.8 cm 
(1.5 in) wide. Although ballistic projectiles were used, the relatively high aspect ratio for these plates could 
possibly allow more global deflection than the specimens impacted in the previously cited studies in this 
section. Results to note are that larger zones of delamination were found to exist for those specimens 
impacted while preloaded. The only exception to this trend was a reduction in the delamination area near 
the back plies on preloaded samples for larger energy impacts. Residual tensile strength was also 
determined for these specimens. It is evident that the tensile preload causes a reduction in residual tensile 
strength by approximately 10 percent for the impact velocities used in this study. 

Data for thinner plates impacted at low-velocity (drop- weight) impacts have been presented 
by Nettles and Lance. 120 Compressive preloads were applied to 16-ply test coupons of T300/934 
and IM7/8551-7 carbon/epoxy with a quasi-isotropic stacking sequence of [0,45,90,-45] S 2- Coupon 
dimensions were 7.6 cm (3 in) wide by 10.2 cm (4 in) gauge length. The preloaded specimens were 
supported over a 6.4-cm- (2.5-in-) diameter circular opening and clamped around its periphery. A design 
of experiments approach was undertaken to try to identity significant variables without testing numerous 
specimens. A total of 1 8 specimens (two replicates of 9) of each material type were tested. Preloads 
ranged from 890 to 44,500 N (200 to 10,000 lb) and impact energies varied between 1 and 16 J 
(0.7 and 1 1.8 ft- lb). In general, it was found that the compressive preloads tended to cause a decrease 
in residual compression strength after impact as the preload was increased, although quantitative values 
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are suspect due to the small number of specimens tested and the variables not being allowed to go to their 
extremes. 

A recent study 121 examined the effect of impact response on biaxially loaded specimens 
of chopped E-glass in a polyester resin matrix. The preloads could be introduced in any combination 
of tension and compression in the two principal axes of the plate. Impacts were performed by a drop- 
weight apparatus at an energy of 21.5 J (15.9 ft-lb) for all specimens. Results measured were indentation 
depth, absorbed impact energy, maximum force of impact, and damage area. The largest level of damage 
was observed to occur for specimens loaded in equal tension-compression (i.e., pure shear). The 
maximum force of impact was largest for the highest values of tension-tension preloading and lowest 
for the largest values of tension-compression (pure shear). A state of biaxial compression produced 
maximum loads of impact between these two extremes. Very little energy was absorbed for specimens 
loaded in tension-tension but as compressive preloads were introduced, the absorbed energy dramatically 
increases. Care must be taken when interpreting these results to continuous fiber-reinforced, layered 
laminates since the damage modes are so very different. 

Analytical treatments of this problem have been undertaken in a few studies. 122-125 In all of these 
studies it was shown that a tensile prestress will decrease the maximum deflection due to impact and will 
reuse the maximum load of impact. Compressive preloads will do the opposite. 
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3. ANALYTICAL 


In this chapter, finite element solutions to the problem of transverse impact of a prestressed plate 
are made. None of these solutions are detailed to the extent that heavy experimental testing need not be 
performed. In keeping with the empirical nature of this work, the analyses presented will be simple and 
produced more to explain observed physical phenomena rather than to predict extent of damage as a 
function of impactor kinetic energy and applied preload. First, an approximate energy method model 
of the problem is presented. A solution of this problem can only be performed numerically and many 
assumptions have to be made, making the finite element technique much more informative and accurate. 
However, some basic physical phenomena can be qualified from the energy method problem before 
it is completely solved. These results can help explain some of the results from the finite element analysis. 

3.1 Energy Method 

The energy involved in the system is assumed to consist of the elastic potential energy stored 
in the deflected laminate and the energy due to the external transverse load. The energy stored in a laminate 
due to bending is given by 
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where b is the width of the laminate, L is the length of the laminate, and w is the transverse displacement 
of the laminate. The D i} are the components of the bending stiffness matrix, x is the cooridinate along the 
length of the plate, and y is the coordinate along the width of the plate. 


The energy stored by midplane (membrane) stretching is given by 


l y-bx-L 
Emem “ ~ / / 

-^-Ox-O 



dxdy , 


( 2 ) 


where the Ay are components of the extensional stiffness matrix. 
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The energy stored due to the externally applied preload, N„ is given by 


I y-bx-L 

ENx = ~ f f 

Z y-0x-0 


where N x is the x-direction stress resultant. 


The energy put into the system by the transverse load, P, is given by 

y-bx-L 4p 

Ep = ~ f / —wdxdy . (4) 

>’=0x-0 Lt> 


The total energy of the system is given by 



Etot m Ebend + Emem + Enx + Ep • (5) 

Assume the deflected surface to be of the form 


w( x, y) m c 0 (x 2 - Lx) 4 sin , (6) 

where C 0 is a constant, satisfies the boundary conditions at x=0 and x=L. The boundary conditions 
at y=0 and y-b restrict the plate from deflecting along these edges when in actuality they are free to move. 
(However, from experimental data to be given later, these edges deflect very little so this approximation 
will suffice for qualitative purposes.) 


Assume the typical values for carbon/epoxy as: 
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(units in pounds and inches) 


and putting equation (6) into equations (1), (2), (3), and (4) and performing the integrations 
(some numerically) with L=16 cm (6.3 in) and y=9.5 cm (3.75 in), then using equation (5) gives 

7. 10 x 10 20 Cq + (8.3 x 10 7 N* + 4.66 x 10 10 )Cq - 1.02 x 10 5 ^Co - Etot , (7) 


with units in pounds and inches. 


From the principle of virtual work, 


dEtot Q 

<?C 0 


( 8 ) 
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Applying this operation to equation (7) gives: 


2.84 x 10 21 C§ + (1.66 x lO 8 ^* + 9.32 x 10 10 )C 0 - 1.02 x 10 5 P . 

At the center of the plate, u^u^, x= U2 and y=b/2. Therefore, 

Co = 1.032 x 10 ^Wftiax > 


(9) 

( 10 ) 


and equation (9) becomes 


P-30 ) 600w^ + (0.168Ax + 94.3)w /na * ■ (H) 

The cubic equation (9) can also be solved for C 0 and then these values placed in equation (6) to determine 
the deflection at any point on the plate. Table 1 gives load (P)-preload (N x ) combinations with the 
associated transverse displacement at the center of the plate (Wim,*). The same values are obtained using 
equation (11). These values are plotted in figure 1 1 for constant values of P with varying N x and in figure 
12 with constant N x and varying P. 

From equation (7), if E to , represents the total impactor kinetic energy put into the system, then if N x 
increases, P increases and C 0 decreases in order to maintain the energy balance of the system Thus, from 
an energy balance formulation of the problem of a transverse load acting on a prestressed plate, the higher 
the prestress, the larger the maximum transverse load and the lower the maximum transverse deflection 
for a given amount of energy put into the system. 
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Prestrain (Microstrain) 

Figure 11. Deflection, u^, versus preload for constant values of transverse load, P from equation (9). 

Load (lb) 


80 160 240 320 400 480 560 



Load (N) 

Figure 12. Load versus displacement for constant values of prestrain from equation (9). 




3.1.1 Formation of Damage 


The above analysis does not take into account any damage that forms in the laminate. From 
equation (5), the energy due to formation of this damage must be included in the total energy. Since 
this energy is lost to the system, it will have a negative value and would tend to decrease the terms 
in equation (11), implying larger deflections will occur for a given transverse load if damage is present, 
thus the “softening” effect. 


3.2 Finite Element Technique 

For this study, a commercially available finite element program called COSMOS/M was utilized. 
The plate was modeled as being statically loaded with a transverse load at the geometric midpoint of the 
plate. Two edges were clamped while the other two were free. Figure 13 shows the plate geometry used 
for the finite element analysis. Large displacement formation was used where the structural stiffness 
matrix was recalculated after each incremental loading of 44.5 N (10 lb). This was necessary to account 
for the changing geometry of the problem since the midplane strains could no longer be ignored; i.e., the 
strain within the laminate midplane is given by: 

du Udw\ 2 

Ex= dx + 2\dx) 


_dv_ Udw \ 2 

£y dy + 2\dy) (12) 

dv du ( dw 
Exy ~ dx + dy + \dx dy) 


where u is displacement in the x-direction and v is displacement in the y-direction. For small displacements 
the final terms in equation (12) can be dropped. However, for displacements greater than approximately 
one-half the plate thickness, these terms can become significant. The finite element large displacement 
formulation performs calculations retaining these terms and recalculating them at each step in the loading 
process. 
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Figure 13. Plate geometry for finite element model. 

The plate was meshed with rectangular composite shell elements 1.016 cm (0.4 in) by 0.914 cm 
(0.36 in) for a total of 160 elements. Each element contained eight layers of orthotropic material with 
properties and directions ([+45, 0,^15, 90] s ) similar to those of the actual specimens to be used. This mesh 
provided adequate aspect ratios for the elements and also produced a grid that was not too course, yet not 
so fine as to extend computational runtime and memory. The transverse load was modeled as a point load 
at the top, centermost node in the mesh. The clamped boundary conditions were achieved by restricting all 
displacements and rotations at all nodes along the sides that were to be clamped. 

Preload was introduced into the laminate by initially displacing one of the clamped ends a certain 
amount in the ^-direction. 

After each run, the desired output response (z-direction displacement, stress in the principal 
material direction of any ply, etc.) could be loaded and plotted. 

3.3 Finite Element Results 


33.1 Transverse Displacement 

The most easily compared result between the finite element (FE) analysis and experimental results 
is the transverse load/displacement relationship. The FE runs were performed for transverse loads of 223; 
445; 668; 890; 1,113; 1,335; 1,558; 1,760; and 2,003 N (50, 100, 150, 200, 250, 300, 350, 400 and 
450 lb) with four prestrains of 0; 2,000; 4,000; and 6,000 he at each transverse load value for a total 
of 36 runs. 

The output from each analysis was presented as a contour plot and the maximum transverse 
displacement (the node under the transverse load) value was recorded. The results from the runs are 
presented in table 2. Figure 14 shows results for a constant transverse load P with varying preload and 
figure 15 for a constant preload N x with a varying transverse load. 
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From the results it can be seen that the maximum deflection decreases with increasing prestrain 
in a nonlinear fashion. The results from the energy balance formulation showed a linear relationship that 
had a much smaller effect (fig. 11). Initially, the added prestrain causes a rapid drop in maximum 
deflection and as the prestrain is increases, this drop becomes less pronounced. As the transverse load 
increases, the absolute drop in maximum deflection is decreased by a very small amount. That is at a 
transverse load of 1,335 N (300 lb), the maximum deflection decreases a total of 2.642 mm (0.104 in) 
as the prestrain increases from 0 to 6,000 fie. For the same increase in prestrain at a transverse load of 
2,003 N (450 lb), the maximum deflection decreases a total of 2.591 mm (0.102 in). 

The load-displacement plots for given prestrain values show that as the prestrain increases, the 
maximum load needed to achieve a given deflection increases by a substantial amount. As the prestrains 
increase, the load-displacement plots become more linear. This was not the case for the energy balance 
results as plotted in figure 2. 

The finite element method predicts a larger influence of the preload on the impact response of the 
laminate than does the energy balance method. Comparisons with experimental results will be made in 
chapter 6. 
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Figure 15. Finite element results for displacement versus load for constant values of prestrain. 
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3.3.2 Shape of Deflected Surface 


From the contour plots of transverse displacements output by the FE program, the shape of the 
deflected surface could be predicted. The surface tended to form a “bowl” shape around the loading point 
rather than deflect as a cylindrical surface (as would be expected for long spans). This precluded any 
closed form analysis of the problem as one of cylindrical bending which would have greatly simplified 
the problem. 

Contour plots for four combinations of transverse load-prestrain values are given in figures 16-19. 
These values represent both no prestrains and high prestrains at a low transverse load and a high 
transverse load. 

As the prestrain increased, the edges along the free sides of the laminate deflected less and even 
began to deflect upward at the highest prestrain, as demonstrated in figures 7 and 9. With no preload, the 
specimen’s free edge did deflect, but at a relatively small amount compared to the maximum deflection, 
as shown in figures 6 and 8. 

3.3.3 Maximum Stresses Within Plies 

The maximum stresses within each ply were also calculated by the finite element method. 

Table 3 gives results (in principal material directions) for the four highest transverse loads for which 
the maximum stresses were calculated. These values are presented to demonstrate general trends 
in the maximum ply stresses with increasing transverse load and prestrain. The maximum stresses 
presented are those at the center of the respective ply. This value was used since most of the stress arises 
from membrane stretching and the values are only an approximation to begin with. 

For the lower plies, the maximum stress increases with both increasing transverse load and 
increasing prestrain, an expected result. However, this trend also holds for the upper plies as well. This 
is to be expected with increasing prestrain but with increasing transverse loads, one might expect the upper 
plies to go more into compression due to the higher bending involved. This demonstrates the significance 
of the large deflection theory as the higher transverse loads add more to the overall membrane-type 
stresses than to the bending stresses. The maximum stresses in the plies in both principal material 
directions decrease from the bottom to the upper plies. Note that only the top two plies ever go into 
compression in the data presented, with these values being very small. Thus, compressive-type failures are 
not expected unless they are due to some mechanism not accounted for in the finite element analysis such 
as contact stresses or dynamic effects. 
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Figure 16. Contour plot of transverse deflection for a transverse load 
of 223 N (50 lb) and no prestrain. 
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Figure 17. Contour plot of transverse deflection for a transverse load 
of 223 N (50 lb) and a prestrain of 6,000 fie. 
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Figure 18. Contour plot of transverse deflection for a transverse load 
of 2,003 N (450 lb) with no prestrain. 
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Figure 19. Contour plot of transverse deflection for a transverse load 
of 2,003 N (450 lb) and a prestrain of 6,000 fie.. 
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Table 3. Finite element results for stresses within plies. 



The strength value in the second principal material direction (-35 Mpa (5,100 psi) as determined 
by 90° tension tests) in the bottom four plies is always exceeded. However, synergistic effects with 
neighboring plies can help keep the ply from failing in this direction at the values obtained by the 90° 
tension tests. It is still expected that much matrix cracking parallel to the fibers in these plies should occur, 
especially in the bottom ply which has the highest stresses and also has only one neighboring ply to help 
“hold together” this ply. 

The breaking stress (-2,700 Mpa (400,000 psi)) in the first principal material direction (along 
the fibers) is never exceeded in any of the plies, thus fiber breakage is not expected to occur for the load- 
prestrain values presented. In fact the stresses calculated are never even one-half that of the breaking stress 
of a ply in the fiber direction. Thus, if fiber breakage is found, it must be due to another mechanism other 
than the stresses set up solely by the transverse load-prestrain combination. 
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4. EXPERIMENTAL 


Before the specifics of this particular research program are given, a few general comments about 
the empirical study of impact to composites are in order; however, direct comparisons between studies are 
very difficult due to a lack of standardization, allowing whatever convenient specimen geometry or sup- 
port system best suits the researcher. However, a hasty rush for a standard can, and would, do more harm 
than good at this point in the evolution of composites, since much remains to be understood about the 
many variables involved in impact testing. Even a standard meant to compare composite materials for 
damage resistance and damage tolerance (such as the Boeing 48-ply CAI standard) can be very mislead- 
ing. For example, in this standard, a 48-ply specimen of a quasi-isotropic stacking sequence, 8 inches 
by 10 inches in size is impacted with 1,500 in- lb of energy per inch of specimen thickness. Assuming 
all impact events are conducted with the same boundary conditions, impactor size, shape, and velocity 
(which they are not), the only response measured is damage tolerance in a structurally compressive test, 
which is inconclusive for assessing a fiber/resin system for impact damage tolerance. It should be noted 
that compression testing a panel of this size is a structural, not material, property test. The most important 
parameter — velocity of the impactor — is not specified, thus high-velocity/low-mass impacts could be used 
to ob tain the required energy level, resulting in very localized damage and a smaller reduction in CAI 
strength than would have been found from a high-mass/low-velocity test. In addition, 48 plies is a 
relatively thick laminate and would not be encountered in many engineering applications, thus bringing 
the issue of scaling into the problem. The thick laminate is called for to minimize the global buckling 
of the plate during compression testing but during the impact event, the flexural stiffness is so high 
that the large interlaminar shear stresses that normally develop in an impact event are minimized. Thus, 
a material with a higher mode II shear strength would not register as being more damage resistant 
than a material with a lower mode II shear strength, although it would for thinner test panels. 

Any standard developed is going to have to examine a range of impact energies. The 1,500 in-lb/in 
called for on the Boeing CAI tests represents just one datum. A material may perforate at this level, giving 
a CAI strength equal to the compression strength of a panel with a hole the size of the impactor drilled 
through it. Another, more damage-resistant/tolerant material may develop a large area of fiber breakage 
(not penetration) associated with delaminations that could produce a much lower CAI strength than the 
less damage-resistant material. At lower impact energies this would have been detected, since perforation 
would not have occurred in the less damage-resistant material and a much larger area of damage would 
have developed whereas the tougher material would probably show a smaller region of damage and 
a higher CAI strength. In other words, it is entirely possible to have material A be more damage- 
resistant/tolerant at one impact energy level than material B, but at a different impact energy level, 
material B can be more damage-resistant/tolerant than material A. 

4.1 Material 

The fiber/resin system chosen for study was IM7 fiber in 977—2 toughened epoxy. This material 
was readily available and represents a very common type of system (high-strength carbon fibers in a 
toughened epoxy) being used on many programs for NASA. The layup configuration chosen for the 
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impact tests was a quasi-isotropic one of [+45,0,^t5,90] s . This configuration is typical of the basic layup 
unit being used for many primary structures on air and spacecraft. The material was processed from unidi- 
rectional prepreg tape as flat panels, using a hot press. The cure cycle consisted of vacuum bagging, then 
heating the panel at 1 .7 °C (3 °F) per minute up to 1 77 °C (350 °F) and holding for 2 hr. Cooling was 
then ramped back down to room temperature at a rate of 2.8 °C (5 °F) per minute. A constant pressure 
of 550 kPa (80 psi) was applied during the entire cycle. The resulting panels had nominal ply thicknesses 
of approximately 0.12 mm (0.0047 in). 

A select series of tests were performed on specimens to determine some key material property 
values. These specimens were manufactured from the same prepreg roll as the impact specimens. The 
cure cycle also remained the same for these specimens. 

4.1.1 0° and 90° Tensile Testing 

Unidirectional laminates of IM7/977-2 were tested to obtain load/strain data so the modulus 
and Poisson’s ratio of the laminae in the fiber direction and perpendicular to the fiber direction could 
be determined. Strain-gauge rosettes were attached to five different specimens of each orientation and 
readings were taken at 445-N (100-lb) increments for the 0° specimens and at 4.45-N (10-lb) increments 
for the 90° specimens. These data were then plotted and a least-squares regression line fit determined the 
measured modulus of elasticity and Poisson’s ratio in the respective direction. The averages were 
En=l,867 Gpa (27.1 msi), v 12 =0.37, E 22 = 9.85 Gpa (1.43 msi), and v 2 i=0.0187. 

4.1.2 ±45° Tensile Testing 

Tensile testing was performed on ±45° specimens to determine the shear modulus of the laminae. 
Twenty-two-ply specimens were tested with strain-gauge rosettes with readings being taken at 223-N 
(50-lb) increments. The stress divided by two versus lei I + \e2\ was plotted and a least-squares regression 
fit on the linear portion of the plot gave an average Gi 2 of approximately 7.03 Gpa (1.02 msi). 

4.1.3 Laminate Testing 

Tensile specimens of the eight-ply [+45, 0,-45,90]$ laminates were also tested with strain-gauge 
rosettes. Readings were taken at 445-N (100-lb) increments and the data were plotted and a least-squares 
fit used to determine the tensile modulus and Poisson’s ratio. The average values were Ex=Ey=60.5 Gpa 
(8.78 msi) and Vi 2 =0.32. 

Strength measurements were also taken on these laminates. It should be noted that edge delamina- 
tion occurred before ultimate fracture. This result is expected as high tensile normal stresses can be set up 
in the region of a free edge on a tensile specimen as outlined by Pipes and Pagano. 126 The ultimate aver- 
age failure strain was 12,875 he which corresponds to an average stress of 779 Mpa (113 ksi). The edges 
began to noticeably delaminate at approximately 689 Mpa (100 ksi). The delamination extended far into 
the specimen width before fracture of the 0° plies. A full analysis and discussion of this phenomenon 
is beyond the scope of this text and the interested reader is referred to other references for a more detailed 
treatment of the subject. 127 ’ 128 
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4.2 Impact Specimens 


The specimens manufactured for impact testing consisted of eight-ply quasi-isotropic samples 
with a stacking sequence of [+45,0,-45,90] s , cut into coupons 9.5 cm (3.75 in) wide and 29.8 cm 
(1 1.75 in) long. The coupons were cut from the cured panel with a diamond-wheel table saw. Fiberglass 
end tabs 1 1.4 cm (4.5 in) wide were bonded onto the ends of the specimen with epoxy adhesive to 
facilitate loading into the pretensioning device. The specimen dimensions are given in figure 20. 
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Figure 20. Specimens used for impact testing. 


4.2.1 Specimen Strain-Gauging 

Each specimen had two strain gauges bonded to it so that a uniform prestrain was ensured. The 
gauges were 1.3 cm (0.5 in) from the specimen’s edges to detect strain gradients across the width (see 
fig. 20). This instrumentation was determined from a series of preliminary tests on samples that contained 
numerous strain gauges on the surfaces of the specimen. From these tests it was found that bending was 
never introduced into the specimen so gauges need only be placed on the same surface. In addition, it was 
found that the strain remained constant along the longitudinal axis of the specimen and the only variations 
in strain were across the width of the specimen. As long as the outermost gauges (1.3 cm (0.5 in) in from 
each edge) had the same strain, the ones between them also had this same strain, thus it was determined 
that only the two gauges nearest the edges and on the same side need be measured to ensure a uniform 
strain field was being set up in the specimen. 

4.3 Pretensioning Device 

The device to apply tensile preloads to the specimens basically consisted of two wedge grips con- 
nected to a load cell and a hydraulic piston. As the piston was pumped with hydraulic fluid it would cause 
the grips to move apart, putting the composite specimen in a state of tension. A sketch of the apparatus 
with the various parts is given in figure 21. 
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The wedge adjustment blocks were first tightened together to lodge the specimen in the grips. 
These blocks also prevented the wedges from moving apart (and stop gripping) during the application 
of the preload. As the preload was applied and the strain gauges were monitored, when one gauge deviated 
from the other, that side reading the lower strain could have its adjustment blocks loosened, allowing the 
specimen to expand slightly along that side, causing the strain to rise on that side. These adjustments were 
made until the final level of strain was reached and the adjustment blocks could be fine-tuned to give 
a nearly uniform prestrain across the width of the specimen. 

4.3.1 Boundary Conditions 

The boundary conditions chosen for this study were clamped-clamped/free-free with the clamped 
edges perpendicular to the applied preload. This geometry was chosen since it was demonstrated in a pre- 
vious program 129 that this condition most closely represented a pressure vessel being impacted. Pressure 
vessels are of vital interest to this study since they are tension-carrying structural members. 


Once the level of prestrain was achieved, the top clamp plates were tightened down with the 
12 bolts, each at a torque of 28 N-m (250 in-lb), to obtain the desired boundary condition. A true clamped 
boundary condition would mean that the hydraulic piston could be released and the specimen would con- 
tinue to carry the prestrain between the region of the clamps. However, in practice this was not the case 
and was never achieved despite trying a number of techniques. The methods that came closest to achieving 
this tended to crush the specimen in the region of clamping; therefore, it was determined that the hydraulic 
piston would be held in place during the impact event. Measurements of the strain before and after impact 
indicated that the specimen did not move significantly in the clamps. The high torque on the clamp-down 
bolts ensured that no bending moments could occur under the clamps. 

4.4 Impact Testing 

Once the desired level of prestrain was achieved, the preloading device with the specimen was 
positioned under a Dynatup 8200 drop tower as shown in figure 22. The height of the crosshead was 
adjusted to give the desired impact energy. In this study, the mass of the crosshead with tup was 2.4 kg 
(5.22 lb). Since the crosshead/guideposts interface does not contain linear bearings or another such 
friction-reducing device, an exact prediction of impact energy from the drop height was impossible. 
Although this friction is present, the measured impact energy is quite precise since the velocity of the 
fallin g crosshead is measured just prior to impact. Once a desired drop height was decided upon, all 
impact tests at that corresponding impact energy level were conducted before the cross head was moved 
again. This ensured uniformity in the impact energy levels used. A total of three impact energies were 
used in this study: 3.4, 4.5, and 6.0 J (2.5, 3.3, and 4.4 ft-lb). These levels induced damage from visible 
to nonvisible. 
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Figure 22. Preload and impact apparatus. 


4.5 Postimpact Testing 

In order to assess the effects of tensile preloads on the damage resistance of the specimens tested, 
three forms of postimpact inspection were used. The first type of inspection was a visual one to see if any 
damage could be detected from either surface of the plate. This was followed by a nondestructive x-ray 
inspection and then a destructive cross-sectional examination of the specimen through the damage zone. 

4.5.1 Visual Inspection 

After each specimen was impacted, it was carefully removed from the preload device and 
inspected for visual damage. No magnifying devices were used; however, if damage could be felt and not 
seen, this was recorded as “visual” damage since no special tools were needed to identify this damage. 
The extent and type of visual damage was recorded for each impacted specimen. The type of damage that 
could be detected consisted of back face (nonimpacted side) matrix splitting, front face (impacted side) 
indentation damage, and small matrix splits being felt on the back surface. 

4.5.2 X-ray Inspection 

For specimens that did contain observable damage on one of the two faces, a small circular dam 
of plumber’s putty was placed around the damage on one of the faces to hold in a zinc iodide (Znl) dye 
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penetrant solution as shown in figure 23. The dye penetrant is opaque to x rays and was allowed to soak 
on the specimen (at a damage site) for at least 24 hr. The solution consisted of 100 g of Znl powder 
in 100 ml of water and 100 ml of Kodak Photoflow ™ which breaks the surface tension of the solution 
to allow full penetration into all cracks and delaminations. 



Figure 23. Dam of plumber’s putty to hold zinc iodide solution over damage area. 

For specimens that did not contain damage on one of the two faces, a small hole of ~0.5 mm 
(0.02 in) was drilled through the specimen at the impact site to allow the dye penetrant to be exposed 
to any through-the-thickness damage that may have occurred within the laminate. 

The specimens were wiped clean of the plumber’s putty and excess penetrant after at least a 24-hr 
soak and exposed to x rays with a Torrex 150D x-ray inspection system manufactured by EG&G Astro- 
physics Research. Polaroid type 55 film was placed directly under the specimen to obtain a permanent 
copy of the x-ray signature of each specimen. The dye penetrant was seen to seep into all matrix delami- 
nations and cracks and was easily detectable on the developed film. 

4.5.3 Cross-Sectional Examination 

Some of the specimens were sectioned through a part of the damage zone in an attempt to gather 
more detailed information about the induced damage. This was particularly helpful in identifying areas 
where fibers were broken, and in identifying the interface of a delamination. 

A technique was developed during this study whereby the location of the cross-sectional cut 
relative to the damage as seen by the x rays could be determined. The x rays, which corresponded to 
the specimens at a 1 : 1 ratio, were photocopied onto film transparency. The clear film now had the image 
of the x ray at a 1:1 ratio with the specimen. Strain gauges that would register on the x rays provided 
accurate guides as to the location that the film needed to be placed to be a surface marker of the internal 
damage. The area of impact on the specimen was painted white to enhance the visibility of the x ray 
on the clear film transparency which was matched to the specimen and applied with double-sided tape. 
The damage zone was then cut from the specimen and two clear acrylic plates were used to sandwich 
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the cut specimen with the x-ray pattern on the surface. These clear plastic plates were bonded to the 
specimen with cyanoacrylate ester glue. The x-ray pattern could be seen clearly through the acrylic plate. 
Any cuts made on the sandwiched specimen could now be associated with the damage, as seen by 
the x-ray inspection technique. 

A microscope equipped with a Polaroid camera was used to examine the cross-sectioned speci- 
mens. Magnifications between x 7.25 and x 160 were possible. 

4.6 Static Indentation Testing 

In order to better understand the load/deflection and in-plane strain response of the transversely 
loaded specimens, a series of static indentation tests were performed using a dial indicator to measure 
transverse displacement and various strain gauges to measure the surface strains at given locations. The 
pretensioning device was placed in an Instron 1125 load frame and the indentor was attached to the upper 
crosshead/load cell. After applying the desired amount of prestrain, the specimen was incrementally 
loaded with the indentor at the center of the specimen (as in the impacts) and transverse displacement 
and strain-gauge readings were taken. This step process was continued until a cracking noise was heard, 
indicating specimen damage. 

The specimens were unloaded in the same incremental fashion to obtain a complete load- 
displacement history of the entire process. This could provide valuable information on how much energy 
is being lost during the process. 

In addition, two static indentation tests were performed in which the specimen was loaded to fail- 
ure. Surface strains were not monitored but load-displacement data were. These two tests were performed 
to gain insight on the load-deflection behavior of specimens with a medium level of prestrain, as well as 
gather information on the static behavior once damage becomes severe in the specimen. 

4.6.1 Static Indentation Strain Gauging 

Strain gauges were placed on the specimens so that information could be obtained about the sur- 
face strains at areas of interest. For practical reasons, strain gauges could not be placed at the center of the 
specimen since the indentor would interfere with it from the top and the dial indicator from the bottom. 

A number of strain-gauge setups were used on the 10 samples tested. Table 4 lists the strain-gauge place- 
ment on each sample along with the rationale behind choosing these locations. 

4.6.2 Loading/Unloading Curves 

As mentioned previously, the transverse load was applied incrementally. After initial readings 
were taken at zero-applied load, the specimen was loaded with 44.5 N (10 lb) and a second reading taken. 
The third reading was taken at 1 1 1.25 N (25 lb) and then readings were taken at 1 1 1.25-N (25-lb) 
increments thereafter until a popping or cracking noise was heard, at which point the specimen was then 
unloaded in 222.5-N (50-lb) increments. The transverse load was applied at a rate of 1.27 mm/min 
(0.05 in/min) between stops. 
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Table 4. Specimens for static indentation tests. 



1. Dimensions in mm. 

2. I denotes Indentor location. 

3. Numbers in parenthesis Indicate strain gauge on opposite side of specimen. 
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Table 4. Specimens for static indentation tests (continued). 



1. Dimensions in mm. 

2. I denotes indentor location. 

3. Numbers in parenthesis Indicate strain gauge on opposite side of specimen. 
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5. RESULTS 


This chapter presents the results of the experimental portion of the work. Since this study is 
empirically based, the material presented is quite detailed, and most of the raw data are presented in the 
appendices. Every effort has been made to make this chapter as complete and thorough as possible. Dis- 
cussion of the particulars of the results and noteworthy comparisons are not made in this chapter but rather 
will be made in chapter 6 which is devoted entirely to a discussion of the results. 

5.1 Instrumented Impact Testing 

The amount of tensile preload on the impacted specimens had a measurable effect on the maxi- 
mum load of impact, duration of impact, maximum deflection, and absorbed energy of the specimen. 

The following sections will present the results from these data. 

5.1.1 Maximum Load of Impact 

As calculated in chapter 3, the higher the tensile prestrain, the higher the impact force should 
be for given impact conditions. This was clearly the case for the specimens examined in this study. 

The results for each of the three impact levels used are shown in figures 24—26. As damage forms 
in the specimen (not accounted for in the analysis in chapter 3), the maximum load of impact tends 
to level off, and even drop some with increasing prestrain (which tends to increase the damage). 



Prestrain (Microstrain) 

Figure 24. Maximum load of impact versus applied prestrain for 3.4 J (2.5 ft-lb) impacts. 
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Prestrain (Microstrain) 

Figure 25. Maximum load of impact versus applied prestrain for 4.5 J (3.3 ft-lb) impacts. 



Prestrain (Microstrain) 

Figure 26. Maximum load of impact versus applied prestrain for 6.0 J (4.4 ft-lb) impacts. 


5.1.2 Duration of Impact 

Figures 27-29 show that as the tensile prestrain increased, the duration of the impact event 
decreased. This variable was not examined in the analysis, but as the maximum deflection decreases 


53 



with increasing prestrain (as will be seen in the next section), it is expected that the duration of impact 
should follow the same trend. Note the nonlinearity in this relationship. 



0 1,000 2,000 3,000 4,000 5,000 6,000 7,000 

Prestrain (Microstrain) 

Figure 27. Total time of impact versus prestrain for 3.4 J (2.5 ft-lb) impacts. 
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Figure 28. Total time of impact versus prestrain for 4.5 J (3.3 ft-lb) impacts. 
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Prestrain (Microstrain) 

Figure 29. Total time of impact versus prestrain for 6.0 J (4.4 ft-lb) impacts. 


5.1.3 Maximum Deflection 


The maximum deflection of impact was predicted to decrease with increasing tensile prestrain as 
determined in chapter 3. Plots of the measured data are shown in figures 30-32. These plots are nonlinear 
in a similar fashion to those from the previous section. The maximum deflection data begins to drop rather 
rapidly with increasing prestrain, then level off at the higher end of prestrains. 
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Figure 30. Maximum center deflection versus prestrain for 3.4 J (2.5 ft-lb) impacts. 
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Prestrain (Microstrain) 


Figure 31. Maximum center deflection versus prestrain for 4.5 J (3.3 ft-lb) impacts. 
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Figure 32. Maximum center deflection versus prestrain for 6.0 J (4.4 ft-lb) impacts. 

5.1.4 Absorbed Energy of Impact 

The amount of energy lost during the impact event is termed “absorbed energy.” The absorbed 
energy versus prestrain plots are given in figures 33-35. For low prestrains the absorbed energy is seen 
to decrease with increasing prestrain up to a critical level, at which point it begins to rise. The smaller the 



56 





impact energy the larger the initial drop in absorbed energy with increasing prestrain. This is thought to be 
a result of compliance in the preloading device. Since a true clamped condition was not practical, at low 
prestrains the specimen could displace a given amount in the fixture, thus losing a given amount of 
energy. As the impact energies increased, this amount of energy became a lower percentage and was less 
noticed. As the prestrain increased, the grips tended to help hold the specimen in a true clamped configu- 
ration and there was less displacement of the specimen, thus less energy lost due to this movement of the 
laminate being tested. 



Prestrain (Microstrain) 

Figure 33. Absorbed energy during impact versus prestrain for 3.4 J (2.5 ft-lb) impacts. 



Prestrain (Microstrain) 

Figure 34. Absorbed energy during impact versus prestrain for 4.5 J (3.3 ft-lb) impacts. 
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Figure 35. Absorbed energy during impact versus prestrain for 6.0 J (4.4 ft-lb) impacts. 

5.1.5 Load-Deflection Curves 

For all of the specimens tested, the load-deflection plots show that a significant amount of energy 
is lost during the impact event, even at the lowest impact energy level used. A typical load/deflection curve 
is shown in figure 36 for the lowest level impact at the smallest prestrain and in figure 37 for the highest 
level impact at the highest prestrain. These two plots represent the two extreme cases of impact tested in 
this study. All of the load/deflection curves are given in appendix A. 
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e 36. Load/displacement curve for lowest level impact at lowest pre 
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Figure 37. Load/displacement curve for highest level impact at highest prestrain. 
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5.2 Visual Examination 


In general, as the preload increased, the amount of visual damage also increased. This corresponds 
to the higher impact force induced into the specimen. The first sign of visible damage was always splitting 
of the matrix between fibers on the back 45° face. Contact stresses produced visible damage at the impact 
site in some of the specimens. Table 5 gives values for all of the specimens tested, along with a descrip- 
tion of the visual damage. 


Table 5. Specimens tested with associated parameters and results. 


Specimen 

No. 

Impact 

Energy 

. UL 

Prestrain 

(me) 

Maximum 
Load (N) 

Maximum 

Deflection 

(nut 

Total 

Time 

(msec) 

Absorbed 
Energy (J) 

Visual 

Damage 

LI 

4.5 

5,182 

1^11 

5.3 

8.06 

3.16 

22-mm split on back face 

L 2 

4.5 

3,229 

1,682 

6.1 

8.70 

2.73 

None 

L3 

4.5 

452 

1,615 

7.6 

10.79 

2.93 

18-mm split on back face 

L4 

4.5 

833 

1,598 

7.6 

10.73 

2.73 

17-mm split on back face 

L5 

4.5 

7,399 

1,664 

4.8 

7.88 

3.81 

1 1 6-mm split on back face 
Dent w/short split, front 

L6 

4.5 

5,564 

1,793 I 

5.3 

8.31 

3.11 

47-mm split on back face 

11/I3A 

4.5 

2,430 

1,700 

6.4 

8.95 

2.52 

17-mm split on back face 

11/138 

4.5 

1,960 

1,660 

6.4 

9.31 

2.66 

18-mm split on back face 

11/14A 

4.5 

7,970 

1,736 

4.8 

7.36 

3.35 

65-mm spilt on back face 
Dentw/small split, front 

11/14B 

4.5 

4,205 

1,825 

5.6 

8.01 

2.51 

None 

11/25A 

4.5 

4,268 

1,829 

5.6 

8.25 

2.92 

30-mm split on back face 

11/27A 

4.5 

6,698 

1,727 

5.1 

8.14 

3.65 

60 -mm split on back face 
Dentw/small split, front 

11/28A 

6.0 

5,384 

2,069 

5.8 

8.13 

4.80 

Split across width 

Dent with short split on front 

11/28B 

6.0 

498 

2,025 

8.4 

10.21 

3.89 

21-mm split on back face 

12/1 2A 

6.0 

7,173 

2,074 

5.3 

8.29 

4.66 

Split across width 
Dent with splits on front 

12/12B 

6.0 

2,132 

2,100 

7.1 

8.86 

3.70 

18-mm split on back face 

12/12C 

6.0 

4,156 

2,123 

6.4 

8.06 

4.26 

45-mm split on back face 
Dent w/small split, front 

12/12D 

6.0 

3,053 

2,096 1 

6.6 

8.59 

4.22 

53-mm split on back face. Dent on front 

12/13A 

6.0 

1,594 

1,967 

7.6 

9.20 

3.70 

20-mm split on back face 

12/138 

6.0 

634 

1,931 

8.4 

10.49 

4.04 

19-mm split on back face 

12/14B 1 

3.4 

6,972 

1,633 

4.6 

7.44 

2.32 

27-mm split on back face 

12/14C 

3.4 

987 

1,335 

6.6 

10.66 

2.14 

None 

12/14D 

3.4 

3,387 

1,477 

5.3 

8.69 

2.07 

None 

! 12/14E 

3.4 

5,012 

1,562 

4.8 

8.09 

2.28 

26-mm split on back face 

12/17A 

3.4 

400 

1,308 

6.9 

1154 

2.24 

None 

12/17B 

3.4 

4,334 

1,531 

5.1 

8.38 

2.09 

19-mm split on back face 

12/17C 

3.4 

6,060 

1,615 

4.6 

7.85 

2.12 

24-mm split on back face 

12/17D 

3.4 

2,045 

1,406 

5.8 

9.50 

1.98 

None 
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5.3 X-ray Inspection 


As expected, the specimens that contained more visible damage also contained more internal dam- 
age, mostly in the form of matrix splitting and delaminations. The specimens that had no externally visible 
damage demonstrated that some delaminating occurred (the small drilled hole allowed the penetrant to 
reach these areas). These delamination zones were relatively small compared to the specimens that had 
visual damage associated with them. A typical x ray of this type is shown in figure 38. Note that matrix 
cracking is associated with the delaminations, a feature seen in all of the specimens. 

Q 


0 0.1 0.2 0.3 0.4 



mm 

Figure 38. Typical x-ray signature for a specimen that exhibits no visible damage. 

For the specimens that did contain externally visible damage, the matrix crack on the back 45° face 
was clearly visible and always had delaminations emanating from the crack. These delaminations tended 
to form a lobe shape and end at a crack in the ply above the 45° ply (the 0° ply). A typical x ray with this 
form of damage is shown in figure 39. 
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As the damage in the specimen became more severe due to higher preloads, damage on the 
impacted side of the specimen would begin to form. For these specimens, the dye penetrant would tend 
to “pool up” and form very dark areas on the x-ray signatures as shown in figure 40. This made detection 
of broken fibers difficult since these breaks are hidden by the large amount of dye penetrant. Cross- 
sectional examination was the method used in this study to detect broken fibers. All of the x-ray signatures 
of damage are shown in appendix B. The static indentation test samples are also included. 
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mm 

Figure 40. Typical x-ray signature for a specimen that exhibits visible damage on the impacted side. 


5.4 Cross-Sectional Examination 

For specimens in which the x rays did not provide sufficient enough detail (usually associated 
with fiber breakage), cross-sectional examinations proved useful. The x rays gave an accurate description 
of the damage formed due to the impact event, but fiber breakage was difficult to detect. When specimens 
were cross-sectioned through the dark areas of the x ray, fiber breakage could sometimes be seen. Cross- 
sectional photographs of selected specimens, along with their x rays showing where the cuts were made, 
are shown in figures 41-49. Note that the 0° fiber breakage areas are always in very close proximity to the 
matrix cracks both above and below this area because these areas are subject to the highest stresses, and 
also, as a matrix crack forms, the component of stress perpendicular to the fibers in that layer is now 
transferred to the layer next to it, in the vicinity of the crack. 
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Figure 49. Cross-section of specimen #12/12C, cut 2. 


Figure 41 is a cross-sectional view of a statically loaded specimen that experienced a maximum 
load of 1,780 N (400 lb). The cut on this specimen was made perpendicular to the 0° fibers in order to 
demonstrate that the long dark areas running in the 0° direction consists of heavy matrix cracking in 
the 0° ply. Note that between delamination zone #3 and matrix cracks #4 and #5, there is much matrix 
damage in the lower 0° ply, which shows up as a dark stripe in the 0° direction on the x-ray signature. 

For the remainder of the cross-sectioned specimens presented, the cross-sectional cut was made 
parallel to the 0° fibers in order to check for fiber breakage in the critical 0° plies. Figure 42 demonstrates 
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that much matrix cracking can exist with no fiber breakage. When fiber breakage is seen, as in figure 43, 
a matrix crack in one of the two adjacent plies is accompanying the 0° fiber breaks. The 0° fiber breaks 
(crack #1) are also seen to be in relatively close proximity to the bottom 445° split (crack #3), that is the 
first and most visible form of damage noted. When the same specimen is sectioned at a point where the 
main +45° and -45° cracks intersect, as shown in figure 44, the bottom 0° fiber breakage is closer to the 
bottom 445° split. For more heavily damaged specimens, fiber breakage in the lower 0° ply is directly 
above the 445° ply split as shown in figure 45. 

Broken 0° fibers in the upper part of the laminate are always seen to be a part of a crack running 
either from the ply above it (the top 445° ply) or below it (the first -45° ply) or both as shown in 
figure 46. The crack in the upper 445° ply (crack #2) and the first 45° crack (crack #1) are joined by 
broken fibers in the 0° ply between them. Figure 47 shows a feature that was seen in many of the 
sectioned specimens. A crack will form from a ply adjacent to the upper 0° ply and transverse part of the 
ply until the crack begins to grow parallel to the fibers, causing the 0° ply to be only partly fractured. This 
crack 

is highlighted by white arrows in the figure. Figure 49 shows a -45° ply crack that starts into the 0° ply 
above it but only breaks fibers a short distance into the 0° ply when the crack is deflected parallel to the 
0° fibers. Note the severe damage in all plies in this cross-sectional view. 

Figure 49 demonstrated the difficulty in trying to identify broken fibers from the x-ray signatures. 
From the cross-sectional view, bottom 0° fibers are clearly broken, yet between cracks #1 and #3 on the 
x-ray signature, there is nothing obvious to indicate that fibers are broken (area #2). This was the case 
for the vast majority of specimens. 


5.5 Static Indentation Tests 

Data from the static indentation tests aided in indicating the shape of the deflected surface as well 
as showing the different strain fields associated with a low versus a high prestrain. The lower preloaded 
specimen demonstrated more nonlinearity at the lower end of loading than the heavily preloaded specimen 
as shown in figures 50 and 51. This is to be expected from the analysis in chapter 3. From the unloading 
data it is seen that a hysteresis is associated with all of the specimens tested at the lower prestrain levels. 
This indicates that there is either more damage being formed in these specimens or more compliance 
in the preloading device at the lower prestrains. Results from the x-ray analysis indicates that more 
damage was not formed, indicating compliance in the fixture at low levels of preload. Complete data 
from the 10 specimens tested are given in appendix C (specimen #2 is not included since the dial indicator 
was not set up properly). 
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Figure 50. Load displacement from a statically loaded specimen at a low pr 
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Figure 51. Load displacement from a statically loaded specimen at a high prestrain. 
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From the in-plane strain measurements as determined from the gauges placed on the specimens, 
as shown in table 4, it was clear that little strain was set up along the edges of the specimen due to trans- 
verse loading, especially on the highly prestrained specimens. Compressive strains due to transverse 
loading were always found in the vicinity of the applied load on the top surface and the strains on the 
opposite side (back face) would be tensile unless the preload was high, in which case they became nega- 
tive. Plots of strain versus applied transverse load for the 10 specimens tested are given in appendix D. 
Note that at regions near the clamped edges (up to 57 mm from a clamped edge for specimens #7 and #8), 
both the top and bottom surface are put into tension due to the transverse load. At high preloads, the 
bottom surface sometimes went into compression from the transverse load (static #5), but this appears 
to be an isolated case. 


5.5.1 Static Indentation Tests to Failure 

Two specimens were prestrained and statically loaded until they could no longer hold a load. 
Only load-deflection data were generated from these tests. The results are shown in figures 52 and 53. 
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Figure 52. Load displacement from a statically loaded specimen (loaded until failure). 
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Figure 53. Load displacement from a statically loaded specimen (loaded until failure). 







6. DISCUSSION OF RESULTS 


This chapter is a discussion of the experimental results obtained in the previous chapter, along with 
results from the analytical analysis performed in chapter 3. 

6.1 Instrumented Impact Testing 


6.1.1 Maximum Load of Impact 

For impact conditions in which the specimen did not sustain heavy damage, the maximum load 
of impact increases linearly with increasing preload. This can be seen for all of the data at the 3.4 J 
(2.5 ft-lb) impact energy level, data up to —5,000 fie prestrain for the 4.5 J (3.3 ft-lb) impact energy level 
and up to —4,000 /re for the 6.0 J (4.4 ft-lb) impact energy level. More scatter exists for the higher two 
impact energy levels due to the development of more severe damage. From equation (7) in chapter 3, it 
was predicted that the maximum load of impact was expected to increase at a linear rate with increasing 
prestrain for a given amount of impact energy. 

6.1.2 Duration of Impact 

The total time of impact is seen to decrease in a nonlinear fashion with increasing tensile prestrain. 
As damage is formed in the specimen, the duration tends to be a little longer due to local “softening” of 
the specimen. Thus, the nonlinearity is due mostly to damage. Note that the total duration of impact seems 
independent of impact energy. This should be the case if the transverse load is approximately proportional 
to the transverse deflection. For a spring-mass system (which can approximate the impact problem), the 
total time of impact (contact of mass on spring) is given by: 


t - 



(13) 


where F is the force exerted by the laminate (spring), x is the displacement, and m is the impactor mass. 
Assuming F <*x (or F=kx, fc=“spring constant”), then equation (1) can be written as: 


■-J? 

which is independent of deflection and therefore independent of impact energy since m and k are fixed 
constants. 


77 



6.1.3 Maximum Deflection 


The maximum deflection of impact decreases in an almost linear fashion with increasing prestrain 
for a given impact level. This linear relationship is not obvious from equation (7). This tends to hold true 
of samples in which incurred damage caused departure from linearity in maximum load and duration data. 
This is due to the “softening” which caused departure from linearity in the maximum load and duration 
data, increasing the deflection, yet the load drops which cancels it out. 

6.1.4 Maximum Absorbed Energy 

These data are seen to decrease with increasing prestrain for the lower values of prestrain, then rise 
with increasing prestrain for the rest of the data. It is predicted that as the preload increases, more energy 
would be lost due to the impact event since higher prestrains tended to cause more damage in the speci- 
men. An explanation for this can be gained from the results of the static indentation tests in which more 
hysterises were formed at the lower prestrains, due to fixture compliance. 

Thus at the low preloads, fixture compliance attributes to energy lost in the system, but as the 
preload increases, this compliance tends to become less severe. At higher prestrain levels, damage begins 
to form in the specimens and the energy lost in the system begins to rise again. 

6.1.5 Load-Deflection Curves 

For higher preloads, the load-displacement plots become more linear up until the point of damage. 
At this point on the curve, the load begins to fluctuate at a high frequency as matrix cracks, delaminations, 
and possibly fiber breakage occurs. 

A large amount of hysteresis is observed for all of the impacted specimens regardless of trans- 
verse load-preload combinations. For specimens that contain significant damage, this result is expected, 
but for specimens with little or no damage, little hysterisis would be expected. This implies that energy 
is being lost to the system in modes other than specimen damage. The particular impact tester used in this 
study does not have linear bearings or other such friction-reducing devices between the falling weight 
crosshead and the guideposts. This has always been suspected of causing the system to lose a fair amount 
of energy in the reaction between the falling crosshead and the guideposts during the impact event. This 
suspicion has never formally been confirmed and documented. 

6.2 Visual Examination 

The first form of noticeable damage that formed was splitting of the matrix along the fibers in the 
bottom +45° ply. This is to be expected from the finite element results, presented in table 3, which show 
these plies develop the highest stresses in the second principal material direction, which are almost twice 
as high as the breaking stress for a lamina. 

The next sign of visual damage was usually associated with contact stresses on the impacted face 
of the specimen. A small dent would become noticeable and as either the prestrain or the transverse impact 
load increased, matrix cracking and fiber breakage was noticed inside the indentation. 
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At relatively high impact levels, broken fibers would begin to protrude from the back face of the 
impacted specimen along a crack in the bottom 445° ply. 

This type of damage progression is typical of fibrous laminates that are allowed to undergo large 
deflections during impact, thus setting up high tensile strains on the back surface of the impacted laminate. 

6.3 X-ray Inspection 

The most notable feature of the x-ray signatures for both the statically loaded and impacted speci- 
mens is that the damage tended to form along the direction of the bottom +45° ply. Delaminations, ema- 
nating from a matrix crack in this ply, were always noted on the x-ray results (unless no damage at all 
formed). This feature is shown schematically in figure 54. The mechanics thought to be behind these 
delaminations are given in the next section. 

6.3.1 Initial Damage Formed 

The initial damage consisted of a matrix crack in the bottom +45° ply with delaminations forming 
as shown in figure 54. On either side of the impact site, the delaminations consisted of two zones, one 
usually smaller than the other. As a matrix crack develops in the bottom +45° ply, a free edge is created 
and the stresses must redistribute themselves so that equilibrium is maintained. This “free edge effect” is 
the same mechanism outlined by Pipes and Pagano. 126 Figure 55 shows schematically how the geometry 
of the problem influences the delamination growth, mainly due to the transverse interlaminar stress, O z , 



Figure 54. Schematic of bottom ply splitting and associated delaminations. 
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Oy 


Face 2 



Face 3 of the ply is carrying the ^-direction stresses, just as it was before the crack was formed. 

The stress distribution is such that the closer a point on the face is to the impact point (I) , the higher the 
jc-direction and y-direction (face 2) stresses are. This is known from the finite element results. The 
moment about a line through point A running in the y-direction must be zero for equilibrium conditions 
to be satisfied. The face 3 stress, Ox, creates a moment about this line. To balance this moment, stresses in 
the ^-direction, Oz, are set up along the crack boundary. Since Ox is largest at the point closest to the 
impact point, I, Oz is also largest at this point. As the crack moves farther away from the line through A 
(i.e., as y increases), the ^-direction stress rapidly decays since Ox gets smaller in this direction and also the 
“moment arm” is getting longer, so to maintain a given moment, the ^-direction force gets smaller. The 
opposite should hold true as the crack gets closer to face 3. The “moment arm” becomes smaller, so the 
stresses need to be larger to maintain a constant moment. However, Ox gets smaller, so the moment 
needed for equilibrium also gets smaller at points further away from the impact point I. 

The above also holds true for Oy which also creates a moment about a line through point B 
running along the ^-direction. This will tend to cause symmetry of Oz along the crack for equal values 
of Ox and Oy. However, since Oy is generally smaller than Ox, especially for heavily preloaded specimens, 
the overall effect is to have slightly larger Oz values near the face 4 end of the crack than at the face 2 end. 
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When face 1 is created, the shear stresses also redistribute themselves to maintain equilibrium. 

The shear stresses on faces 2 and 3, Txy, are no longer balanced by shear stresses on face 1, since it is free, 
or on face 4, since it can be considered infinitely small. New interlaminar shear stresses T zy and Tzx are 
created to balance the moments due to the shear stresses. Since Txy is largest nearest the point of impact 
(as it was for Ox and Oy), Tzx and Tzy will also be largest nearest the impact point, I. 

6.4 Cross-Sectional Examination 

Large matrix cracks and delaminations could easily be detected by x-ray analysis but fiber 
breakage was difficult to detect except by cross-sectional examination. Cross-sectional examination 
showed that fiber breakage could occur in a laminate at a point that was not obvious from the x-ray 
signature. A large amount of matrix damage was always associated with neighboring plies of a ply with 
broken fibers. These areas of high crack density tended to cause the dye penetrant to “pool-up” and create 
dark areas on the x-ray signatures which would obscure the broken fiber damage. 

The cross-sectional analysis also revealed that fiber breakage in the bottom 0° ply differed from 
that in the upper 0° ply. The fiber breakage in the bottom 0° ply always ran completely through the ply 
but fiber breakage in the upper 0° ply did not. In many instances the fiber fracture would only travel one- 
half or one-third the thickness of the lamina and then deflect parallel to the fibers, running through the 
matrix in the upper 0° ply. The bottom 0° ply fiber breakage was always associated with large delamina- 
tions between this ply and the bottom +45° ply. The fiber breakage in the upper 0° ply was characterized 
by the propagation of the crack containing the fiber breakage into one, or if the ply was completely frac- 
tured, both adjacent plies through the matrix material in these adjacent plies. Delamination was sometimes, 
but not necessarily, associated with these cracks. 

The failure modes of the two 0° plies can be expected to differ, based on the maximum ply stress 
values given in table 3. The upper 0° plies are subjected to small stress values that are not even within an 
order of magnitude of the fiber-breaking stress for these plies. Thus, the fiber breakage must come from 
another stress mechanism setup, other than that by preload or membrane and bending stresses arising 
from the transverse load. Since dents were observed on the impacted surface of many of these laminates, 
it is assumed that contact stresses are responsible for fiber breakage in the upper 0° ply. The stresses 
predicted for the lower 0° plies are much larger (although they never exceed one-half the breaking stress) 
and with added stress from cracked plies transferring stress to the main load-bearing 0° fibers, it is likely 
that these fibers simply break in tension (causing the complete fracture of the ply). 

6.5 Static Indentation Tests 

The static indentation tests showed that a tensile preload increases the flexural stiffness of the 
laminate, especially at lower levels of transverse load, at which point the induced membrane stresses due 
to large deflection are not as large. 
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6.5.1 Comparison With Energy Balance and Finite Element Models 

Figure 56 shows load-deflection curves for specimens prestrained at 4,000 [xe as predicted 
by equation (9) and the finite element analysis. Data from a static indentation test are included for 
comparison. It is clear that the finite element method comes closer to predicting the actual load- 
displacement behavior of the laminate. 

Figure 57 shows prestrain-deflection data as calculated by the two models, along with data 
from the static indentation tests at a constant transverse load of 223 N (50 lb). The finite element method 
matches up extremely well with the experimental data. The curve as predicted by the energy balance model 
(equation (9)), greatly underestimates the effect of prestrain on deflection and also predicts a linear rela- 
tionship which does not agree with experimental data. 

Displacement (mm) 



Displacement (in) 

Figure 56. Load-deflection data comparison between equation (9), the finite element method, and static 

indentation data for a prestrain of 4,000 fxe. 
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Figure 57. Prestrain-deflection data for a constant transverse load of 223 N (50 lb). 

The data from the in-plane strain measurements agreed qualitatively with the finite element results 
and showed that at high preloads, little deflection or strain is induced into the laminate near the free edges 
due to the transverse load. 

A result that was found experimentally that was not predicted was the compressive strains that 
appeared on the top surface of the laminates due to the transverse load during the static indentation tests. 
From table 3 it was predicted that as the transverse load increased, the stress (or strain) on the top surface 
would become more positive due to the large membrane stresses set up in the laminate. One cause for this 
could be the previously mentioned contact stresses which are not accounted for in the model. 

6.5.2 Comparisons of Impact and Static Indentation Curves 

From all accounts, the impacts in this study were of a quasi-static nature with no contributions 
from material inertia or material vibrations. The conditions for a quasi-static impact event, low target 
mass, and a large impactor mass were certainly met in this study. Thus, it is expected that the static 
indentation data should be very similar to the impact data. 

Figures 58-60 show load-deflection data for three different levels of preload for both the static 
indentation and impact tests. These data are seen to match fairly well, especially at the lower end 
of transverse deflection and at lower preloads. The static indentation data show a larger load for a given 
deflection at the higher end of deflection. Note that the data from the impact tests show a small drop 
in stiffness along the load-deflection curve whereas this drop is not seen in the data from the static 
indentation tests. In addition, figure 59, which has static indentation data for a specimen with heavy 
damage, shows that the static indentation specimens can sustain a slightly higher load before the first 
sudden drop in load is observed. In figure 60, no damage has occurred on the static indentation specimen 
at 1,780 N (400 lb), but obvious damage has occurred in the impacted specimen. 
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Figure 58. Comparison of load-deflection curves for a static indentation specimen 
and an impact specimen both preloaded at -228 N/cm (130 lb/in). 
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Figure 59. Comparison of load-deflection curves for a static indentation specimen 
and an impact specimen both preloaded at -1,300 N/cm (742 lb/in). 
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Figure 60. Comparison of load-deflection curves for a static indentation specimen 
and an impact specimen both preloaded at -2,600 N/cm (1,500 Ib/in). 

It was noticed during the static indentation tests that when a specimen was at a load that could 
cause damage formation, damage did not always immediately occur in the specimen. As the load was held 
constant and a displacement reading was being taken, the specimen would sometimes produce a loud pop 
which indicated the formation of damage after a time interval at a given load. 

From x-ray results on specimens that experienced similar maximum transverse load-preload 
combinations, some impacted specimens showed more damage than the specimens that were statically 
loaded, although the differences were not great. Most of these specimens contained very similar x-ray 
signatures. 

Thus some subtle differences can be noted for impacted and statically loaded specimens in this 

study. 
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7. CONCLUSIONS AND RECOMMENDATIONS 


7.1 Overview 

This chapter contains conclusions that can be drawn from this study and recommendations 
for future studies of relevance. 


7.2 Conclusions 

The following are the major conclusions of this study: 

• For all other impact parameters which are the same, a tensile prestress tends to increase the maximum 
force of impact and the amount of damage formed. A decrease is seen in the duration of impact and in 
the maximum transverse deflection due to impact. 

• The types of damage formed are matrix cracking (splitting), delaminations, and fiber breakage. The 
first sign of damage is matrix splitting in the bottom +45° ply. Large delaminations emanate from 
the matrix split and tend to run between this split to a split in a neighboring ply. Smaller delaminations 
form directly under the impact zone due to the high contact stresses set up in this area. 

• Large deflection plate theory must be taken into account for analysis in which the plate deflects more 
than one-half its thickness. The induced midplane stretching from the large deflections rapidly 
dominate the load/displacement response of the plate. 

• Finite element techniques do a good job of predicting the elastic response of a transversely loaded 
plate, but the introduction of damage greatly complicates the problem. 

• X-ray inspection can easily detect matrix cracks and delaminations; however, fiber breakage 
is difficult to detect. 

• Cross-sectional examination revealed that fiber breakage did occur, even though the finite element 
analysis predicted that the maximum stresses along the fibers in a ply would not even be one-half the 
breaking strength of the ply. 

• Fiber breakage in the lower 0° ply was dominated by tensile stresses whereas fiber breakage 
in the upper 0° ply was dominated by contact stresses and showed a complex failure mode. 

• Statically loaded specimens and impact loaded specimens show the same general behavior 

if the target mass is small and the impactor mass high, but subtle differences do exist, especially at 
high tensile prestrain values. 
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7.3 Recommendations 


The following recommendations are for related programs that may be conducted in the future 

concerning impact damage of composite laminates: 

• The gray area between “quasi-static” and dynamic cases must be accounted for since many impacts 
may fall into this regime. 

• Finite element analysis must be handled with care once damage is incorporated into the model. 

• The inclusion of the x-displacement values, u, in the energy balance solution should be incorporated 
for more accuracy involving the influence of the preload. 

• The effects of compressive preloads should be studied, since by the arguments given in section 6.3.1, 
if a crack forms, oz should be compressive and help to suppress delamination growth. 

• Nondestructive evaluation (NDE) techniques need to be examined extensively since, in the field, this 
will ultimately be the tool that determines the amount of damage present due to an impact event. 
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8. SUMMARY 


Tensile preloads tend to increase the maximum force of impact and decrease the duration 
of the impact event and the maximum transverse deflection for low- velocity impacts. Higher tensile 
preloads also tend to cause more damage to form in the specimen due to the higher impact forces and 
normal tensile stresses, oz, that are set up along a back face split that helps to form a delamination. 

Finite element analysis can predict the material response well but once damage is induced into 
the specimen, the problem becomes much more complex. 

Fiber breakage is difficult to detect from x-ray signatures of the damage zone due to the high 
matrix crack density associated with such a specimen, causing the dye penetrant to “pool-up” and hide 
the fiber breaks. 

Contact stresses are very important, even when the specimen is flexible and allowed to deform 
many times its thickness. Static indentation and “quasi-static” impact tests show the same basic behavior 
and damage formation, but subtle differences do exists. 


88 



APPENDIX A— LOAD-DEFLECTION DATA FROM INSTRUMENTED IMPACT TESTS 
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APPENDIX B— X-RAY SIGNATURES OF STATIC INDENTATION 
AND IMPACTED SPECIMENS 
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Static indentation specimen # 4 




Static indentation specimen # 7. 



Static indentation specimen # 9. 
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Static indentation specimen # 10. 
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APPENDIX C— LOAD-DISPLACEMENT CURVES FROM STATIC INDENTATION 

TESTS 
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APPENDIX D— STRAIN GAUGE READINGS FROM STATIC INDENTATION TESTS 
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